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ABSTRACT 


An  array  of  l86  sheet-metal  collector  pads  was  placed 
in  a  concentric  circular  pattern  about  ground  zero  (GZ)  of  a 
near- surface  nuclear  explosion.  The  detonation,  with  a 
yield  of  approximately  0.5  kt  (4.53  X  10^  kg),  occurfed  at  a 
depth  of  58  cm  in  desert  alluvium.  After  the  shot,  when 
residual  radioactivity  permitted,  samples  of  ejecta  anrt  dust 
(throwout)  deposited  on  these  collectors  were  obtained, 
weighed,  and  analyzed  for  particle  size.  One  hundred  sam¬ 
ples  were  thus  recovered  at  distances  ranging  from  50  to 
600  meters  from  GZ.  In  ballistic  trajectories,  large  parti¬ 
cles  were  found  to  have  traveled  the  greatest  distances; 
deposition  of  very  fine  particles  on  the  outermost  rings  is 
postulated  to  have  occurred  at  least  partially  as  a  result 
of  the  base  surge.  Mass  (weight)  distribution  of  throwout 
per  Tinit  area  is  tabulated  and  shown  graphically  as  a  func¬ 
tion  of  radial  distance  from  GZ,  and  comparisons  with 
similar  experiments  are  made. 

Effects  of  yield,  depth  of  charge  burial,  and  other 
variables  upon  the  ejecta  equation  are  examined,  as  are 
means  of  expressing  the  eqviation  itself.  Methods  of 
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predicting  both  ejecta  areal  density  and  thickness  are  pre¬ 
sented  and  discussed,  as  are  relations  between  certain 
mechanics  of  crater  formation,  crater  volimie,  and  ejecta 
volume . 
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dv,  dt,  etc. 
D  ,  D 


NOTATION 


a  coefficient  for  the  dimensionless  expression 
of  the  ejecta  eqn/ition 

a  coefficient  foj-  the  ejecta  equation  written 

2 

for  dimensionless  distance  from  GZ,  kg/m 

particle  diameter 

apparent  crater  depth,  meters 

effective  diameter  of  an  ejecta  particle 

depth  of  true  crater,  meters 

first  derivatives  of  the  quantities  shown 

drag  force 


function 


gravitational  force 


2 

acceleration  due  to  gravity,  n/sec 
height  above  ground  surface,  meters 


unit  vector  in  x-direction 


unit  vector  in  y-direction 

a  coefficient  in  the  hyperbolic  expression  of 
the  ejecta  eq'uation,  6  =  k 


a  constant 


unit  vector  in  z-direction 
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an  exponent  in  the  ejecta  eqtiation  denoting 


decay  of  ejecta  hreal  density  with  increasing 

radiaJ.  distance  from  GZ 

mass  of  ejecta  pkrticle 

effective  mass  of  ejecta  particle 

mass  of  the  appak:‘ent  crater  volume,  kg 

mass  of  ejecta,  kg 

an  exponent 

common  ratio  of  a  geometric  progression 
apparent  crater  tadius,  meters 
true  crater  radios,  meters 
radius  vector  describing  the  position  of  a 
particle  of  ejedta;  i.e.  the  magnitude  and 
direction  of  itd  slant  range  from  the  origin 
radial  distance  (range)  from  GZ,  meters 
radial  distance  from  GZ  to  the  crater  lip 
crest,  meters 

the  horizontal  distance  traveled  by  an  ejecta 

particle  from  its  original  position  in  the 

crater,  meters 

time  after  detonation,  sec 

thickness  of  ejdcta  deposition,  meters 
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u 


V 

V 


V 

o 


V 

s 


V  . 

XI 

V 

X  ,  y 


V 


the  velocity  of  an  ejecta  particle  with 
respect  to  the  ground 
velocity  at  any  given  tine 
terminal  velocity  for  a  particle  in  inter¬ 
mediate  flow  conditions 
initial  velocity 

terminal  velocity  for  a  particle  in  streamline 
flow  conditions 
terminal  velocity 

terminal  velocity  of  a  particle  in  turbulent 

flow  conditions 

initial  horizontal  velocity 

velocity  at  any  time  parallel  to  the  x-  or 

y-axis 

the  velocity  of  an  ejecta  particle  with 
respect  to  the  air  (or  medium)  mass 
apparent  crater  volume,  m 

that  part  of  the  volume  of  a  crater  caused  by 
compression  and  plastic  flow,  m 

3 

volume  of  ejecta,  m 

volume  of  ejecta  at  seme  value  of  R/r  ,  m'^ 
true  crater  volume,  m 
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wind  velocity 


charge  weight  or  yield  expressed  as  eqtiiva- 
lent  weight  of  TNT,  kt  or  kg 
axes  of  a  three-dimensional  coordinate  system 
scaled  position  of  charge  with  reference  to 
the  earth-air  interface 

scaled  depth  of  burst ,  ft/kt^^^  ’  ^  cr  n^/kg^'^^ 
the  specific  value  of  ^  which  describes  the 
direction  of  a  constant  horizontal  wind, 
degrees 

specific  weight  (density)  of  ejecta,  kg/m 
ejecta  areal  density;  mass  (weight)  distribu¬ 
tion,  kg/m^ 
a  coefficient 

angle  of  elevation  of  ejected  particle,  meas¬ 
ured  from  the  vertical,  degrees 
angle  of  elevation  of  ejected  particle  at  the 
instant  of  ejection,  degrees 
angle  between  an  ejecta  particle  velocity 
vector  and  the  horizontal 

the  angle  between  a  particle  velocity  vector 
and  the  horizontal  at  the  instant  of  ejection 
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K 


K. 

1 


K 

S 


K 


t 


X 

It 

P 

P 

( 

a 


a  constant  associated  vrith  the  dimensionaJ. 
expression  for  drag  force  on  a  particle 
following  a  trajectory  in  a  fluid 
the  constant  <  as  applied  to  intermediate 
flow  conditions 

the  constant  K  as  applied  to  streamline  flow 
conditions 

the  constant  K  as  applied  to  ttirbulent  flow 
conditions 

mean  free  path  of  air  molecules 
viscosity  of  a  given  medium 
a  constant,  3 *1^16 
particle  density 
density  of  the  air 

strength  properties  of  the  cratered  medium 
angle  between  the  z-axis  and  the  xz  plane, 
degrees 


oo  infinity 

One  dot  (•)  denotes  the  derivative  with  respect  to  time 

✓  N 

.  Two  dots  (••)  denote  the  second  derivative  with 
respect  to  time. 
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CHAPTER  1 
INTRODUCTION 

1.1  OBJECTIVE 

The  objective  of  Project  1.5  was  to  determine  the  dis¬ 
tribution  of  crater  ejecta  (throwout)  resulting  from  the 
detonation  of  a  low-yield,  near- surface  nuclear  device  in 
desert  alluvium,  including  the  distribution  by  weight 
(or  mass)  and  grain  size. 

1.2  BACKGROUND 

The  mass  movement  of  earth  materials  by  surface  and 
near- s\arf ace  nuclear  detonations  is  of  military  significance 
primarily  because  of  its  effects  upon  nearby  military  in¬ 
stallations,  especially  those  "hardened"  to  survive  near- 
misses  of  nuclear  warheads.  Damage  caused  by  ejecta  strik¬ 
ing  or  covering  antennas,  missile  site  components,  entrances 
to  protective  structures,  etc.,  presents  a  definite  hazard. 
The  behavior  of  ejected  material  is  also  of  interest  in 
peacetime  projects  involving  nuclear  energy,  such  as  the 
excavation  of  canals,  harbors,  etc. 

Tne  history  of  ejecta  experiments  shows  a  wide  variety 
of  environmental  conditions.  Because  of  the  variation  in 


CONFIDENTIAL 

RESTRICTED  DATA 


media,  data  which  are  ultimately  to  be  employed  to  support 
scientific  hypotheses  will  probably  req\xire  treatment  by 
statistical  methods,  establishing  probabilities,  randomness 


or  variation,  interdependence  of  variables,  mean  values, 
standard  deviations,  and  similEU*  statisticsil  parameters. 
There  is  also  a  need  for  a  more  quantitative  definition  of 
the  parameters  governing  ejecta  deposition,  as  well  as  their 
association  with  scaled  depth  of  burst,  crater  radius  and 
depth,  etc.  The  solution  of  these  problems  will  require  a 
large  amoiint  of  data.  The  Johnie  Boy  shot  provided  a  means 
of  atigmenting  the  information  available  from  the  studies 
discussed  below  and  was  of  particular  significance  because 
of  its  near- sxarf ace  geometry. 

Prior  to  Johnie  Boy,  43  cratering  shots  had  been  made 
which  furnished  usable  crater-ejecta  data.  Several  studies 
have  been  made  (References  1  through  lO)  of  the  areal  depo¬ 
sition  and  particle  size  distribution  of  throwout  from  cra¬ 
ters  produced  by  explosive  charges  ranging  in  size  from  a 
few  hundred  pounds  to  0.5  kt  of  high  explosive  (HE)  and  from 
0.42  to  100  kt  of  nuclear  explosive  (NE).  These  studies 
have  in  some  cases  been  included  as  primary  projects  and  in 
other  cases  as  secondary  investigations.  Charge  positions 
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have  varied  from  that  of  a  contact  btirst  (resting  t^jon 
ground  surface)  to  a  scaled  depth  of  0.68  io/'aq  ' 

(216  ft/kt^^^)  in  several  types  of  soil  and  rock.  Forty- 
shots  have  been  HE  and  three  NE.  Synopses  of  these  in¬ 
vestigations  are  provided  in  the  following  paxagraphs . 

1.2.1  Engineering  Research  Associates,  Inc.  (ERA). 

ERA  dust  (ejecta)  studies,  conducted  in  Utah,  included  media 
of  sandstone,  limestone,  dry  sand,  dry  clay,  and  -wet  clay, 
and  were  performed  with  trinitrotoluene  (TNT)  charges  vary¬ 
ing  in  size  from  1.45  X  10^  to  1.45  X  10^  kg  (320  to  320,000 

1/^ 

pounds)  placed  at  scaled  depths  of  0.l4  to  0.40  nv/kg  ' 

(44  to  127  ft/kt^^^).  In  the  sand  and  clay  studies,  samples 
were  obtained  at  distances  from  ground  zero  (GZ)  of  125  to 
500  feet  on  smaller  shots  and  to  distances  of  approximately 
3,100  feet  on  the  larger  shots.  In  the  sandstone  tests,  the 
sampling  range  was  50  to  1,000  feet  for  the  small  shots  and 
750  to  5,000  feet  for  the  largest.  The  single  limestone 
test  utilized  sampling  gages  50  to  250  feet  from  GZ.  Re¬ 
sults  of  these  studies  are  summarized  in  Reference  1,  and 
additional  tab-ulation  and  plotting  of  the  data  are  included 
in  Reference  10. 
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Two 


1.2.2  Suffield  Experimental  Station  (SES). 

contact  shots,  one  of  5  tons  and  the  other  of  100  tons  of 
TNT,  were  fired  on  the  Watching  Hill  site  at  the  Suffield 
Experimental  Station  in  Ralston,  Alberta,  Canada,  and  throw- 
out  studies  were  made  (References  3  and  4).  The  surface  de¬ 
posits  at  this  site  consist  of  silt  and  clay,  underlain  by 
soft  strata  of  clay,  shaTe,  sandstone,  and  gravel  to  a  depth 
of  300  feet. 

1.2.3  Nevada  Test  Site  (NTS)  Studies— HE  Shots.  Eight 
HE  cratering-eJecta  experiments  have  been  conducted  at  the 
U.  S.  Atomic  Energy  Commission's  Nevada  Test  Site:  Scooter, 
Stagecoach  I  and  II,  White  Tribe,  Jangle  HE-2,  and  Buckboard 
11,  12,  and  13.  Except  for  the  Buckboard  shots,  which  were 
in  basalt,  and  White  Tribe,  which  was  in  caliche,  the  deto¬ 
nations  were  conducted  in  desert  alluvium.  All  were  btiried 
shots  except  White  Tribe,  which  dealt  only  with  contact 
bursts.  Jangle  HE-2  was,  however,  quite  near  the  surface. 
References  2,  5?  6,  8,  and  9  report  these  events. 

1.2.4  NTS  Studies — NE  Shots.  In  addition  to  Johnie 
Boy,  three  NE  cratering-ejecta  experiments  have  been 
conducted  to  date  at  NTS:  Teapot  Ess,  Danny  Boy,  and 
Sedan  (References  2,  7,  and  lO).  These  three  shots  in¬ 
volved  deeply  buried  charges  and  therefore  are  not 
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directly  comparable  to  Johnie  Boy. 

1.2.5  Simmary  of  E.jecta  Data.  Appendix  A  siDimarizes 
in  tabular  form  all  usable  cratering-ej'ecta  experiments 
known  to  the  authors;  the  tabiilation  includes  also  the  re- 
sialts  obtained  from  Johnie  Boy.  Figures  1.1  and  1.2  are 
selected,  normalized  ejecta  curves  from  Reference  10. 

1.3  THEORY 

1.3»1  General.  As  shown  in  Reference  10,  there  is  a 
definite  q\ialitative  relation  between  ejecta  distribution 
and  other  cratering  parameters,  such  as  scaled  depth  of 
bvirst,  crater  radius,  depth,  volume,  etc.  Thus^  ejecta 
parameters  are  siffected  by  many  of  the  same  factors  that 
affect  true  and  apparent  crater  formation,  as  well  as  by 
weather  (wind  and,  to  a  lesser  degree,  precipitation). 

That  is,  apparent  crater  volume  is  dependent  upon  the  size 
of  the  true  crater  and  upon  the  volume  of  the  fallback  ma¬ 
terial.  The  apparent  crater  results  from  (l)  ejection  of 
material  from  the  crater,  (2)  flowage  of  material  into  the 
upthrust  region  which  lies  beneath  the  crater  lip,  and  (3) 
compression  of  the  in  situ  material  that  lies  just  beyond 
the  botindary  of  the  true  crater. 
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Charge  parameters  include  size,  energy  density,  shape, 
velocity  of  detonation,  the  pressure-time  history  of  the  de¬ 
veloped  gas  bubble,  thermal  properties,  efficiency,  depth  of 
burial,  and  effective  or  scaled  depth.  Medium  parameters 
include  density,  strengths,  homogeneity,  compressibility, 
stress  wave  attenuation  properties,  viscosity,  and  directly 
related  gravity  effects. 

Particle  size  distribution  of  the  ejecta  and  its  rela¬ 
tion  to  air  drag  amd  wind  effects  are  important  elements  in 
determining  the  distance  particles  of  a  given  size  will 
travel  from  the  crater.  For  large  particles  (missiles), 
ballistic  trajectories  are  determined  by  the  initial  veloc¬ 
ities  imparted  to  the  particles,  the  angles  of  elevation  of 
the  trajectories,  and  the  drag  resistance  of  the  atmosphere. 
For  small  particles,  the  settling  (terminal)  velocities  in 
relation  to  wind  become  important.  Thus  the  velocity  field 
generated  by  the  explosion  in  the  cratered  material  and  the 
degree  of  comminution  largely  determine  the  particle  trajec¬ 
tories  and  the  distances  the  particles  will  travel  from 
their  initial  positions  within  the  crater. 

1.3.2  Cratering  by  Shallow  Bursts.  The  interrelation 
between  cratering  mechanisms  varies  with  scaled  depth  of 
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burst  and  the  characteristics  of  the  cratered  medium.  For 


shallow  biursts  in  sand  and  fine-grained  soil,  the  mechanism 
is  dominantly  that  of  viscous  flow,  and  it  is  assumed  to  be 
similar  for  tuff  and  alluvium.  Reference  11  reports  a  com¬ 
puter  program  developed  for  a  theoretical  model  of  a  2-Mt 
surface  burst,  assuming  that  the  input  required  (l)  the  his¬ 
tory  of  the  nuclear  explosion,  (2)  an  equation  of  state  for 
the  earth  medium,  and  (3)  an  estimation  of  the  energy  par¬ 
titioning.  A  hydrodynamic  model  was  e]T5)loyed  in  the  early 
stages,  and  gravity  was  not  a  consideration  in  the  crater 
formation.  The  usual  hydrodynamic  equations  for  shock  phe¬ 
nomena  were  employed  to  describe  the  conservation  of  mass, 
momentum,  and  energy,  with  a  pressure  function  of  two  vari¬ 
ables,  density  and  specific  internal  energy.  Effects  of 
heat  conduction  eind  viscosity  were  neglected,  although  it 
has  been  postulated  (Reference  12)  that  viscous  effects  and 
a  consequent  l/2-power  scaling  law  may  dominate  the  forma¬ 
tion  of  craters  that  are  produced  by  very  large  surface 
bursts.  The  shock  energy  was  governed  by  Hugoniot  condi¬ 
tions.  Results  obtained  on  the  two-dimensional  model  thus 
established  were  found  to  be  in  reasonable  agreement  with 
experimental  observations. 
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Although  the  model  was  developed  for  a  2-Mt  burst  in 
ttiff,  the  results  should  apply  in  principle  to  the  Johnie 
Boy  experiment.  Theoretically,  the  pressures  and  tempera¬ 
tures  of  a  small  burst  should  be  the  same  as  those  of  a 
large  burst,  anaJLagous  to  the  process  of  detonation  of  chem¬ 
ical  explosives  (HE),  where  the  peak  temperature  and  pres¬ 
sure  are  independent  of  the  charge  size,  provided  optimum 
detonation  velocity  is  achieved.  It  would  also  be  expected 
that  the  vectors  in  the  velocity  field  would  be  geometri¬ 
cally  similar  for  identical  scaled  depths  of  btirst  of  the 
same  explosive,  with  adjustments  made  for  the  actual  and 
gravity- induced  differences  in  the  properties  of  the  media. 

As  a  result  of  the  very  high  presstires  generated  by  a 
nuclear  explosion,  allwium  not  only  flows  but  is  consider¬ 
ably  coK^jressed.  It  is  estimated  in  Reference  2  that  com¬ 
pression  of  the  surrounding  material  accounted  for  4l  per¬ 
cent  of  the  apparent  volume  of  the  Sedan  crater.  However, 
Reference  13  indicates  that  the  size  of  an  explosion  crater 
is  not  materially  affected  by  volumetric  changes  of  the  soil 
within  which  the  crater  is  formed.  It  appears  logical  to 
assume  that  the  soil  is  temporarily  compressed,  but  recovers 
most  of  its  original  vol\mie,  and  that  the  apparent  crater 
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volume  comprises  principally  contributions  of  ejecta  volume 
and  plastic  flowage  of  the  soil. 

The  intense  pressure  of  the  detonation  undoubtedly 
causes  both  fusion  and  pulverization  of  s'urrounding  soil  and 
rock.  Small  particles  of  ejecta  are  more  severely  affected 
by  air  resistance,  and  settling  velocities  relative  to  wind 
velocity  determine  the  distances  to  which  small  particles 
will  be  carried.  While  no  exact  quantitative  relations  can 
be  derived  theoretically  concerning  ejecta  distribution, 
velocity  fields,  and  wind  effects,  a  discussion  of  the  basic 
principles  involved  is  given  in  Appendix  B  in  order  to  offer 
a  partial,  qualitative  explanation  and  to  indicate  the  areas 
of  need  for  more  ftmdajnental  research. 
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Figure  1.1  Comparison  of  ERA  and  Buckboard  ejecta  curves 
(Reference  lO) . 
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Figrire  1.2  Comparison  of  Stagecoach,  Scooter, 
and  Danny  Boy  ejecta  curves  (Reference  lO) . 
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CHAPTER  2 


PROCEDURE 

2.1  SITE  LOCATION 

The  Johnie  Boy  event  was  conducted  in  area  l8  of  the 
NTS  (see  Fig\ire  2.1). 

2.2  ENVIRONMENTAL  CONDITIONS 

The  test  area  was  a  relatively  flat  portion  of  an  allu¬ 
vial  fan.  The  soil  was  gravelly  sand,  containing  numerous 
cohhles  and  exhibiting  some  cementation;  there  was  little 
variation  in  its  composition  to  a  depth  of  at  least  80  feet, 
the  deepest  boring  made  by  personnel  of  the  Soils  Division, 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES),  dur¬ 
ing  the  soils  survey.  Soil  density,  measured  near  GZ,  was 
approximately  1,870  kg/m^  (ll7  Ib/ft^).  The  rather  sparse 
vegetation  in  the  area  consisted  mostly  of  sagebrush. 

The  measured  properties  of  the  soil  (Reference  l4)  at 
the  Johnie  Boy  site  include  stratification,  soil  type,  water 
content,  mineral  content,  density,  and  particle  size  distri¬ 
bution.  Samples  taJcen  in  the  vicinity  o'f  GZ  showed  gravelly 
sand  with  cobbles  to  a  depth  of  30  feet,  with  sand  and  occa¬ 
sional  cobbles  to  the  full  depth  of  the  hole.  Depth  of  the 
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stratum  with  largest  percentage  of  cobbles  varied,  being 
generally  greatest  near  GZ.  Because  of  difficulties  caused 
by  cobbles  during  the  drilling  of  the  upper  20  feet  of  the 
soil,  samples  from  this  region  are  not  considered  tr\jly 
representative . 

Mineralogical  analysis  showed  that  the  larger  sieve 
sizes  were  composed  of  dense,  light  gray  to  pink  volcanic 
rocks,  vesicular  porphyries,  and  glassy  particles.  Almost 
all  of  the  particles  were  coated  with  a  calcareous  covering. 
In  general  the  samples  consisted  principally  of  rocks  rang¬ 
ing  from  soda  dacite  to  rhyolite,  with  textures  ranging  from 
porphyries  to  tuffs  and  highly  glassy  rocks  with  flow  struc¬ 
tures  shown  in  elongated  pebbles. 

The  particle  size  curves  in  Reference  l4  show  much  more 
fine  material  below  30  feet,  with  an  increase  of  plus  No.  10 
material  above  this  level.  For  the  near- surface  detonation 
of  the  Johnie  Boy  device,  it  would  be  expected  that  most  of 
the  ejecta  would  come  from  above  the  30-foot  depth. 

2.3  EXPERIMENTAL  ARRAY 

2.3.1  Geometry.  A  circular  pattern  of  throwout  col¬ 
lectors  (see  Figure  2.2)  was  chosen  due  to  (l)  the  somewhat 
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tmpredictable  wind  spectrxnn,  and  (2)  the  random  circumferen- 
tiaJ.  variation  in  throwoTit  deposition  known  to  have  occvirred 
in  previous  experiments.  The  greatest  radial  distance  (600 


meters)  was  chosen  to  represent  30  predicted  crater  radii, 
the  maximum  distance  at  which  recovery  of  identifiable 
ejecta  samples  was  anticipated.  The  minimum  distemce  (50 
meters)  was  selected  to  fail  slightly  outside  the  predicted 
crater  lip.  The  choice  of  seven  circular  rings  of  collector 
pads  was  an  arbitrary  one,  based  largely  upon  limitations  of 
time  and  funds.  With  the  foregoing  conditions  fixed,  it  was 
assumed  that  throwout  deposition  would  decrease  with  in¬ 
creasing  radial  distsuice  from  GZ  as  a  geometric  progression 
of  common  ratio,  r  .  Therefore,  r^  =  ,  and  r  1.5  , 

the  mxiltiplier  used  to  determine  the  radial  distances  of  the 
rings  between  50  and  600  meters.  The  number  of  collectors 
on  each  ring  was  approximately  proportional  to  the  circum¬ 


ference  of  the  ring,  begin.iing  with  a  minira\im  of  12  on  the 
three  inner  rings.  A  total  of  I86  collectors  was  included 
in  the  array. 


2.3.2  Placement  and  Numbering  of  Collectors.  Each 


collector  consisted  of  a  2.50-by-3.00-foot  (0.70  m  )  sheet- 


metal  pad.  Stirvey  for  the  collector  layout  was  provided  by 
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Holmes  and  Narver,  on-site  archi  tect-engineers,  who  marke^ 
all  collector  stations  by  means  i  f  coded  stsl^es.  Where  nec 
essary,  a  small  area  was  cleared  at  each  station,  and  col¬ 
lectors  were  roughly  leveled  and  oriented.  In  some  cases, 
as  when  stations  fell  on  or  near  reads,  the  collector  posi¬ 
tion  was  moved  to  one  side  fax  enoi:,^h  to  prevent  its  being 
disturbed.  Spikes  driven  through  predrilled  holes  Jn  +.he 
metal  pads  served  to  secure  collectors  in  place.  In  addi¬ 
tion,  collectors  on  the  three  innermost  rings,  which  were 
subjected  to  the  greatest  blast  and  shock,  were  set  in  con¬ 
crete.  The  rings  were  lettered  in  alphabetical  order,  be¬ 
ginning  with  A  as  the  closest  to  GZ,  and  collectors  on  each 
ring  were  numbered  in  clockwise  fashion,  number  1  being  the 
first  station  east  of  true  north.  Figure  2.3  is  an  aerial 
view  of  the  Johnie  Boy  site  during  preshot  construction,  in 
which  a  portion  of  the  Project  1.5  layout  is  visible.  To 
gxiard  against  preshot  acciunulation  of  debris,  regular  in¬ 
spections  of  the  collector  rings  were  made  until  just  prior 
to  shot  time. 

2.4  TEST  CONDITIONS 

The  Johnie  Boy  shot  was  detonated  at  0945  hours. 
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2.5  RECOVERY  OF  TISOtfOUT  SAKPLES 

It  v«ji  Jjqportant  that  the  particulate  deposition  on  the 
collectors  be  recovered  as  soon  as  possible,  in  order  to 
■inlmize  the  effects  of  weather  on  the  fine  siaterials.  Ini¬ 
tially,  residual  radiation  ssade  it  necessary  that  exasiina- 
tion  of  the  collectors  be  lisiltcd  to  the  outer  rings  and  to 
the  southern  portion  of  the  array.  Recovery  began  on  D  +  1 
(12  July)  and  continued  intermittently  (as  reuiiation  levels 
permitted)  until  26  Jxily,  at  which  tine  135  collectors  had 
been  examined  and  73  samples  recovered.  The  remainder  of 
the  recovery  work  was  accomplished  on  11  October.  Of  the 
l86  collectors,  100  samples  were  obtained,  21  collectors 
were  destroyed  or  buried  by  debris  so  as  to  make  recovery 
impracticable,  and  3  collectors  were  found  to  have  been  dis- 
t\irbed  so  as  to  render  any  samples  on  them  useless.  The 
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rcsiainder  (62)  contained  no  identifiable  deposition.  Sev¬ 
eral  of  the  close-in  pads  had  been  dislodged  by  the  force  of 
the  blast,  but  the  concrete  bases  were  found  in  place,  and 
samples  were  taken  from  these.  The  initial  effor?:  ’..'s.s  to¬ 
ward  recovery  of  the  maximum  number  of  specimens  in  the 
shortest  time  possible.  Therefore,  samples  were  sealed  in 
metal  containers,  which  were  marked  and  left  in  place.  They 
were  later  removed  to  a  shelter  about  4,000  feet  from  GZ, 
where  initial  weighing  was  accomplished  and  data  sheets  were 
prepared.  The  samples  were  then  resealed  and  stored  for 
later  sieve  auialysis.  Figures  2.5  and  2.6  show  the  recovery 
and  weighing  operations.  The  extent  of  sample  recovery  is 
graphically  illustrated  in  Figure  2.7* 

2  .6  SIEVE  ANMiYSES 

WES  Soils  Division  personnel  performed  sieve  analyses 
on  99  throwout  samples  (one  was  lost  in  handling)  in  Feb¬ 
ruary  1963*  For  this  purpose,  the  sealed  samples  were 
transported  to  CP-2  at  NTS,  the  location  of  Reynolds  Elec¬ 
trical  and  Engineering  Company's  Rad-Safe  office.  Here  they 
were  stored  and  processed  in  a  shelter  and  trailer  set  aside 
for  that  purpose  and  under  conditions  designed  to  facilitate 
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radioactivity  control  and  saitty  measures.  U.  S.  standard 
sieves  and  a  mechanical  shsikel  were  used.  A  second,  more 
accurate  weight  was  also  ob^i^ied  on  each  specimen,  iinall 

t 

samples  were  analyzed  in  their  entirety,  while  standard 
quartering  procediires'were  employed  to  reduce  large  samples 
to  a  more  convenient  size. 

1 ' 

2.7  INVESTIGATION  OF  (^RATER  LIP 

In  June  19^3 ,  in  conjunction  with  a  study  of  crater 
meas\irements,  excavations  were  made  in  the  region  of  the 
crater  lip  to  define  the  interface  between  the  upthrust 
original  gromd  surface  and  the  ejecta  overburden.  Figure 

2.8  shows  the  area  of  this  investigation,  which  was  taken 
to  a  distance  of  46  meters  (150  feet)  from  GZ.  At  the 
greater  distances,  small  pits  were  dug,  and  the  ejecta 
thicknesses  were  measured.  Where  the  crater  lip  was  well 
defined,  trenches  were  dug  extending  north,  east,  and  south 
of  the  crater  (to  the  west,  excavation  had  already  occurred 
on  another  project),  and  the  original  grovmd  surface  was 
mapped  by  means  of  conventional  survey.  The  profile  of  the 
crater  lip,  including  the  ejeCta  overburden,  was  also  sur¬ 
veyed.  Apparent  sliding  and  folding  of  bedding  planes  near 
the  crater  edge  made  identification  of  the  ground  surface- 
ejecta  interface  uncertain  in  seme  places.  Observations 
made  by  WES  Geology  Branch  personnel  are  included  as 
Appendix  C. 
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Figxire  2-3  Aerial  view  of  Johnie  Boy  site,  showing  part  of 
Project  1.5  layout.  Activity  is  preshot  survey,  construc¬ 
tion,  etc.  Note  concrete  crosses  used  for  ground  reference 


stations.  Apparent  asymmetry 
tortion  in  photograph. 


in  layout  is  caused  by  dis- 
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NOTE:  YIELD  =  0.5  KT. 

Fig\zre  2.4  Shot  geometry,  Johnie  Boy. 
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Figure  2.6  Initial  weighing  aind  recording  of  samples.  The 
Howe  scale  had  a  300-ponnd  capacity  in  l/ll-pound  gradtia- 
tions;  the  Soiltest  scale  (in  use)  had  a  35-po'und  capacity 
in  0.01-pound  graduations. 
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Close-in  investigation  of  ejecta  deposition 
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CHAPTER  3 


RESULTS 

3 . 1  MASS  DISTRIBUTION 

Figure  3-1  is  an  aerial  stereopair  of  the  Johnie  Boy- 
crater,  showing  the  lip  and  the  general  pattern  of  ejecta 
distribution.  Weights  of  recovered  samples  are  given  in 
Table  3-1.  Recovery  of  ejecta  samples  from  the  A  ring  was 
poor  due  to  the  effects  of  the  blast  and  to  the  ejecta 
thickness,  which  made  it  difficult  to  locate  the  collector 
pads.  The  recovery  for  rings  B  through  G  was  satisfactory; 
however,  some  of  the  collectors  on  F  and  most  of  those  on  G 
had  no  readily  discernible  ejecta  deposited  on  them.  The 
sample  weight  data  are  also  presented  graphically  on  an 
ejecta  distribution  map  (Figure  3 "2).  Although  at  shot  time 
there  was  a  7-mph  wind  from  195  degrees  true  azimuth,  it  had 
no  apparent  effect  upon  ejecta  distribution. 

The  thickness  of  ejecta  deposition  near  the  crater  lip 
is  shown,  together  with  the  upthrust  of  the  ground  surface, 
in  Figure  3*3-  Tables  3-2  and  3.3  give  ejecta  thicknesses 
in  the  region  of  the  crater  lip.  These  data  were  obtained 
on  three  radials  only,  and  they  provide  a  general 
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representation  of  ejecta  profiles,  which  actvially  exhibited 
a  considerable  circumferential  variation.  For  detailed  pro¬ 
files  of  the  lip  crest,  see  Appendix  C.  Preshot  and  post¬ 
shot  contoiir  maps  were  employed  to  calculate  the  total  ver¬ 
tical  change  in  siorface  elevation  near  the  crater  (see 
Figvire  3*^)-  The  contours  in  Figure  3-^  reveal  the  amount 
of  upthrust  and  throwout;  they  were  obtained  by  referring 
all  elevations  to  a  horizontal  plane  through  the  (preshot) 
GZ. 


3.2  PARTICLE  SIZE  DISTRIBUTION 

Each  sample  was  screened  through  U.  S.  standard  screen 
sizes  1- 1/2-inch,  3/4- inch.  No.  10,  and  No.  200,  or  with 
Nos.  8  and  16  in  place  of  No.  10  (see  Table  3-4).  For  pur¬ 
poses  of  analysis  and  graphical  comparison,  the  screen  frac¬ 
tions  obtained  with  Nos.  8  and  16  were  converted  to  an  ap¬ 
proximate  distribution  corresponding  to  that  for  No.  10  only 
by  placing  one-quarter  of  the  No.  16  fraction  with  the  No.  8 
fraction  and  the  remaining  three-quarters  with  plus  No.  200. 
The  adjusted  No.  8  fraction  was  then  assumed  to  define  the 
No.  10  fraction. 

The  cumulative  distribution  for  each  ring  was 
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calcxilated  using  the  total  weight  of  each  screen  fraction 
for  all  samples  in  the  ring.  For  the  F  and  G  rings,  how¬ 
ever,  samples  F-1  and  G-32  were  omitted  from  the  calcula¬ 
tions  because  their  weights  were  very  large  compared  with 
the  other  samples  in  those  rings,  and  their  inclusion  would 
have  biased  the  values  of  the  cumulative  distributions. 

Their  omission  is  permissible  because  the  total  cumulative 
distribution  gives  weight  in  proportion  to  the  weights  of 
the  individual  screen  fractions,  and  any  statistical  treat¬ 
ment  to  determine  significauit  samples  requires  that  the  sam¬ 
ples  obtained  be  representative  of  the  population  from  which 
they  are  drawn.  As  each  ring  was  selected  to  represent  a 
stratimi  which  was  presumably  homogeneous  in  character,  those 
samples  which  do  not  fall  into  the  pattern  of  homogeneity 
may  be  rejected.  Although  there  are  mathematical  tests  to 
determine  the  acceptability  of  sauries,  the  general  charac¬ 
ter  of  the  data  here  does  not  seem  to  justify  the  applica¬ 
tion  of  detailed,  rigorous  mathematical  criteria. 

The  cumulative  particle  size  distributions  of  the 
composite  ring  samples  and  those  for  the  in  situ  alluvium 
(from  Reference  l4)  are  shown  in  Figure  3- 5- 
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3 . 3  EJECTA  PHOTOGRAPHY 

Six  photographs  of  the  fireball  and  early  trajectory 
history  of  the  Johnie  Boy  shot,  taken  for  a  project  con¬ 
ducted  by  Edgerton,  Germeshausen,  and  Grier,  Inc.  (EG&O), 
show  the  phenomena  accompanying  the  near- surface  burst  at 
times  ranging  from  0.48  to  5*62  sec  after  detonation 
(Figure  3-6). 
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TABLE  3.1  MASS  DISTRIBUTION  OF  EJECTA 


Station 

Radial 
Distance 
from  OZ 

Field 
Weight  of 
Specimen 

Metric 
Wel^t  of 
Specimen 

Ejecta 

Distribution 
per  Unit  Area 

Station 

Radial 
Distance 
from  GZ 

Field 
Weight  of 
Specimen 

Metric 
Weight  of 
Specimen 

Ejecta 
Distribution 
per  Unit  Area 

meters 

lb 

ke 

meters 

lb 

kg 

kg/m*" 

A- 3 

50 

12.80 

5.82 

8.35 

E-18 

265 

0.06 

0.03 

0.04(-) 

A-5 

28.00 

12.73 

18.26 

E-19 

0.22 

0.10 

0.14 

E-20 

0.14 

0.06 

0.09(-) 

B-1 

75 

35-5 

16.14 

23.16 

E-21 

0.10 

0.05 

0,07 

B-5 

2.70 

1.23 

1.76 

E-22 

0.07 

0.03 

0.04(+) 

B-6 

9.50 

4.32 

6.20 

B-9 

22.50 

10.23 

14.68 

E-23 

265 

0.21 

0.10 

0.l4 

B-10 

12.00 

5.45 

7.82 

E-24 

0.77 

0.35 

0.50 

E-25 

0.73 

0.33 

0.47 

C-1 

115 

7.10 

3.23 

4.63 

E-26 

0.18 

0.08 

0.12 

C-3 

3-30 

1.50 

2.15 

E-27 

0.10 

0.05 

0.07 

C-5 

2.90 

1.32 

1.89 

C-6 

1.12 

0.50 

0.72 

E-26 

265 

1.81 

0.82 

1.18 

C-8 

0.55 

0.25 

0.36 

E-29 

0.51 

0.23 

0.33 

E-30 

0.48 

0.22 

0.32 

C-9 

115 

4.10 

1.86 

2.67 

C-10 

79.0 

35.90 

51.51. 

F-1 

400 

4.10 

1.86 

2.67 

c-u. 

4.61 

2.09 

3.00 

F-2 

0.04 

0.02 

0.03 

F-3 

0.09 

0.04 

0.06 

D-1 

175 

4.30 

1.95 

2.80 

F-U 

0.01 

0,01 

O.Ol(-) 

D-2 

14.90 

6.77 

9.71 

F-6 

0,10 

0.05 

0.07 

D-3 

5.60 

2.54 

3-65 

D-U 

3.22 

1.46 

2.10 

F-7 

400 

0.23 

0.11 

0.16 

D-5 

1.99 

0,90 

1.29 

F.8 

0.04 

0.02 

0.03(.) 

F-10 

0.07 

0.03 

0.04(+) 

D-6 

175 

o.4o 

0.18 

0.26 

P-U 

0.02 

0.01 

0.01(+) 

D-7 

0.4l 

0.19 

0.27 

P-12 

0.10 

0.05 

0.07 

D.8 

0.24 

0.11 

0.16 

D-9 

5.23 

2.37 

3.40 

F-13 

400 

0.13 

0.06 

0.09(.) 

D-10 

6.08 

2.74 

3.93 

P.l4 

0.09 

0.04 

0.06 

F-15 

0.06 

0.03 

0.04{.) 

D-U 

175 

2.74 

1.24 

1.78 

P-16 

0.06 

0.03 

0.04(+) 

D.12 

0.31 

0.14 

0.20 

P-18 

0.05 

0.02 

0.03 

D-14 

0.25 

0.11 

0.16 

D-15 

1.28 

0.58 

0.83 

P-19 

400 

0.02 

0.01 

0.01(+) 

D-16 

0.62 

0.28 

0.40 

F-20 

0.04 

0.02 

0.03(-) 

P-21 

0.05 

0.02 

0.03 

D-17 

175 

9.81 

4.45 

6.39 

F-22 

0.04 

0.02 

0.03 

D-18 

2.68 

1.22 

1.75 

F-23 

0.06 

0.03 

O.C4(.) 

D-19 

9.00 

4.09 

5.87 

D-20 

8.35 

3.80 

5.45 

F-25 

400 

0.02 

0.01 

0.01(+) 

P-27 

0.08 

0.04 

0.06 

E-1 

265 

3-95 

1.80 

2.58 

F-28 

0.13 

0.06 

0.09(-) 

£-2 

0.17 

0.08 

0.11 

P.29 

0.03 

0.01 

0.01(+) 

E-3 

1.10 

0.50 

0.72 

F-30 

0.05 

0.02 

0.03 

e-4 

1.92 

0.87 

1.25 

E-7 

0.78 

0.35 

0.50 

P-32 

400 

0.92 

0.42 

0.60 

F-33 

0.07 

0.03 

0.04(+) 

e-8 

265 

0,12 

0.05 

0.07 

F-35 

0.03 

0.01 

0.01(+) 

E-9 

0,02 

0.01 

O.Ol(-) 

F.36 

0.06 

0.03 

0.04(+) 

£-10 

0.3)* 

0.15 

0.22 

F-39 

o.o4 

0.02 

0.03 

E-U 

0.05 

0.02 

0.03 

F-40 

0.08 

o.o4 

0.0& 

E-12 

3.94 

1.79 

2.57 

G-32 

600 

0.13 

0.06 

o.09(-) 

E-13 

265 

0,09 

0.04 

o.o6(-) 

G-37 

0.15 

0.07 

0.10 

E-14 

0.02 

0.01 

0.01(+) 

G-40 

0.01 

0.01 

O.Ol(-) 

E-15 

0.06 

0.03 

0.04(+) 

G-42 

0.02 

0.01 

0.01(+) 

E-16 

0.13 

0.06 

0.09(-) 

G-44 

0.02 

0.01 

0.01(+) 

E-17 

0.08 

0.04 

o.o6(-) 

G-55 

0.02 

0.01 

0.01(+) 

G-59 

0.02 

0,01 

O.Ol(-) 
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TABLE  3.2  EJECTA  THICKNESSES  AT  SELECTED  STATIONS 


True  Azimuth 

Ejecta  Thickness 

27*43  meters 
(90  feet) 
from  GZ 

36*57  meters 
(120  feet) 
from  GZ 

45.72  meters 
(150  feet) 
from  GZ 

degrees 

meters 

meters 

meters 

22.5 

0.U6 

0.06 

0.09 

45.0 

0.21 

0.15 

0.15 

67.5 

0.18 

0.27 

0.12  ' 

90.0 

0.55- 

0.55 

0.15 

112.5 

0.30^ 

0.21 

0.12 

135.0 

o.ko 

0.12® 

0.09 

157.5 

0.30 

O.UO 

0.40 

180.0 

0.79 

0.61 

0.15 

202.5 

0.3*^ 

0.12 

0.06 

225.0 

0.18 

0.18 

0.09 

2U7.5 

0.21 

0.09 

0.09 

270.0 

0.30 

0.24 

0.18 

292.5 

0.12 

b 

0.18 

315.0 

0.3*^ 

0.12 

0.15 

337.5 

0.61 

0.24 

0.30 

360.0 

0.49 

0.30 

0.12 

^  Area  disturbed;  thickness  estimated. 
Area  disturbed;  no  meas\irement  taken. 
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TABLE  3.3  EJECTA  THICKNESS  THROUGH  CRATER  LIP 


North  (0  Degree 

True  Azimuth  frcan  GZ) 

East  (90  Degree 

Tn-e  Azimuth  from  GZ) 

South  (180  Degree 

True  Azimuth  from  GZ) 

Distance 
frcm  GZ 

Ejecta 

Thickness 

Distance 
from  GZ 

Ejecta 

Thickness 

Distance 
from  GZ 

Ejecta 

Thickness 

meters 

meters 

meters 

meters 

meters 

meters 

23-77 

1.83 

21.95 

0.30 

22.16 

0.30 

24.38 

1.92 

22.86 

0.34 

22.86 

1.04 

25.91 

1.34 

24.38 

0.37 

24.38 

1.49 

27.43 

0.49 

25-91 

0.40 

27.43 

0.79 

28.96 

0.64 

27.43 

0.55 

30.48 

0.34 

30.48 

0.49 

30.48 

0.61 

36.57 

0.61 

33-52 

0.49 

32.00 

0.61 

45.72 

0.15 

36.57 

0.30 

33-52 

0.52 

38.10 

0.37 

36.57 

0.55 

39-62 

0.15 

39.62 

0.15 

45.72 

0.12 

45.72 

0.15 
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denotes  adjustaent  to  include 


TABI£  3-**  (Concluded) 
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Figure  3-1  Aerial  stereopair  of  the  Johnie  Boy  crater, 
showing  the  lip  and  general  pattern  of  ejecta  distribution. 
For  three-dimensional  effect,  use  standard  stereoscope. 
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Figixre  3-2  Ejecta  distribution  contours. 
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Figure  3-^  Postshot  ejecta  surface  elevations  reduced  to  a 
horizontal  reference  plane. 
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Figure  3.6  Sequence  of  selected  EG&3  photographs  of.Johnie 
Boy,  showing  fireball  and  ejecta  at  designated  times  after 
detonation. 
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CHAPTER  4 


DISCUSSION 

4.1  RELIABILITY  OF  DATA 

4.1.1  General.  Within  the  scope  of  its  planned  ob- 
Jectives,  the  Johnie  Boy  Project  1.5  provided  adeqtiate  in¬ 
formation  for  an  accwate  representation  of  the  ejecta  dis¬ 
tribution  except  perhaps  for  large  particle  sizes.  For  the 
purposes  of  engineering  analysis,  and  considering  budgetary 
limitations,  saniple  collectors  were  placed  in  maximimi  num¬ 
bers  and  sizes.  A  layout  designed  to  yield  sufficient  data 
for  rigorous  statistical  and  mathematical  treatment  would 
have  increased  costs  considerably.  The  Johnie  Boy  ejecta 
study  utilized  l86  sampling  stations,  as  compared,  for  exam¬ 
ple,  with  272  for  the  Sedan  st^ldy.  Thus,  in  relation  to  the 
size  (yield)  of  the  test  shot,  the  Johnie  Boy  throwout  study 
represents  one  of  the  most  detailed  projects  of  its  kind 
conducted  to  date. 

The  close-in  samples  could  not,  because  of  objection¬ 
able  radiation  levels,  be  collected  until  several  weeks 
after  the  event,  and  the  amount  of  fine  material  in  each  of 
these  samples  was  unquestionably  affected  by  wind  and 
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precipitation;  however,  the  amount  deposited  or  removed  by 
winds  was  considered  to  be  negligible  percentagewise.  No 
attempt  was  made  to  control  the  moisture  content  of  the  sam¬ 
ples  during  processing;  however,  in  the  arid  climate  of 
Nevada,  this  is  probably  of  minor  importance.  A  few  samples 
were  wet  when  recovered,  but  these  dried  out  before  final 
processing  was  begun. 

The  question  of  identification  of  ejecta  should  also  be 
considered.  Near  the  crater  edge,  overturning  and  apparent 
sliding  of  bedding  planes  accounted  for  most  of  the  lip 
height.  Overtiirned  planes  could  be  readily  identified, 
principally  by  their  inverted  root  systems,  and  were  treated 
as  ejecta.  Sliding  of  the  planes  was  more  difficult  to  de¬ 
tect,  and  these  sections  were  not  considered  ejecta,  al¬ 
though  it  is  recognized  that  they  contributed  to  the  crater 
volume.  Farther  from  GZ,  there  were  found  on  some  collec¬ 
tors  particles  too  large  to  be  attributed  to  postshot  depo¬ 
sition  by  wind,  and  yet,  by  their  low  radiation  levels,  ap¬ 
parently  not  from  the  crater.  It  seems  likely  that  air 
turbulence  accompanying  the  blast  caused  these  depositions 
by  merely  blowing  or  rolling  particles  along  the  ground  from 
one  point  to  another.  There  seemed  to  be  no  practical  way 
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of  separating  particles  according  to  their  origin,  so  ni  i 
were  treated  as  ejecta. 

4.1.2  Statistical  Considerations.  For  possible 
theoretical/statistical  treatment,  the  deposited  ejecta  mass 
may  be  considered  as  a  population  composed  of  particles  of 
varying  sizes  at  increasing  distances  from  GZ.  The  statis- 
td.cal  parameters  are  weight  per  unit  area,  distance  from  GZ, 
and  particle  size  distribution.  Closely  related  factors  are 
the  particle  size  distribution  of  the  in  situ  alluvitm  and 
the  drag  effects  of  the  air  on  particle  trajectories.  Even 
when  it  is  assumed  that  the  detonation  produces  considerable 
fragmenting  of  in  situ  particles,  the  esctent  of  which  would 
be  difficult  to  evaluate,  the  in  situ  distribution  is  be¬ 
lieved  to  be  closely  related  to  that  of  the  ejecta. 

Altho\agh  the  collector  array  for  Johnie  Boy  was  symmet¬ 
rical  with  respect  to  GZ,  for  statistical  purposes  it  may  be 
assumed  that  a  random  sainpling  of  ejecta  was  obtained.  In 
the  planning  of  experiments  of  this  type,  there  is  no  way  to 
predict  the  circumferential  variation  of  ejecta  deposition, 
and  a  symmetrical  collector  array  is  therefore  best  siiited 
to  obtain  a  random,  representative  group  of  samples.  In  or¬ 
der  that  these  samples  be  equally  representative,  it  is 
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desirable  that  those  on  a  given  ring  be  of  approximately  the 
same  weight.  This  is  particularly  necessary  where  an  arith¬ 
metical  average  is  involved.  As  mentioned  in  Section  3*2, 
two  saii5)les  were  not  included  in  the  averages  for  outer 
rings  because  of  their  excessive  weights. 

4.2  PARTICLE  SIZE  DISTRIBUTION 

Particle  size  distributions  by  ring  (Figure  4.l)  depict 
graphically  the  sorting  that  results  from  trajectory  mechan¬ 
ics  in  air.  The  amounts  of  all  large  material  (No.  10  and 
above)  except  +l-l/2-inch  increased  in  proportion  to  the 
distance  traveled.  Although  the  larger  size  particles  rep¬ 
resent  a  fairly  large  portion  of  the  in  situ  alluvium,  their 
percentage  distribution  by  number  is  small.  It  can  readily 
be  shown  that  the  area  covered  by  the  collectors  in  the 
outer  rings  was  so  small  compared  to  the  few  large  particles 
propelled  to  these  distances  that  the  probability  of  one  of 
them  landing  upon  a  collector  was  remote.  Truly  representa¬ 
tive  sampling  would  have  required  collectors  with  total 
areas  varying  in  proportion  to  distances  from  GZ;  thus, 

tailing  the  collectors  on  ring  A  as  standard,  the  G  ring 

2 

would  have  required  a  sampling  area  of  more  than  100  m  , 
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or  l44  individual  collectors.  Using  this  criterion  for  the 
entire  sampling  array  would  have  required  approximately  400 
collectors.  Thus  a  modified  procedure  based  on  the  above 
criterion  was  employed  by  which  a  convenient  collector  size 
was  made  the  standard,  and  the  numbers  of  collectors  were 
varied  so  as  to  provide  a  sampling  array  that  was  economi¬ 
cally  manageable  within  the  funds  provided  for  the  study. 

If  the  number  of  particles  per  unit  weight  of  material 
rather  than  weight  only  is  considered,  together  with  the 
fact  that  a  one -cubic -inch  particle  is  equal  to  1,000 
l/lOth-inch  particles  and  is  also  equal  to  4.03  X  10'  par¬ 
ticles  of  minus  No.  200  mesh  size,  it  is  readily  seen  why 
the  weight  percentage  of  larger  particles  may  be  much  too 
low  for  samples  obtained  at  the  greater  ranges  where  small, 
uniform-size  collectors  are  employed.  Thus,  while  the 
+3/4- inch  material  increases  out  to  400  meters,  there  was 
(except  for  sample  G-32)  none  of  this  size  e^iecta  found  on 
the  collectors  at  600  meters.  It  can  be  reasoned,  there¬ 
fore,  that  sampling  such  as  used  in  this  experiment  did  not 
give  a  true  measure  of  the  amount  of  coarse  material  which 
was  actually  ejected  to  the  distances  at  which  the  collec¬ 
tors  were  placed.  The  sorting  indicated  did,  however,  give 
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a  qualitative  indication  of  the  effects  of  air  drag  in  de¬ 
termining  the  distances  to  which  ejecta  will  travel,  both 
where  the  trajectories  are  completely  ballistic  and  where 
the  settling  occurs  at  terminal  velocities. 

■'’articles  in  the  size  range  less  then  2  mm  will  quickly 
reach  terminal  velocities  lower  than  1  w/sec  and  would  be 
expected  to  fall  near  the  crater  if  there  were  no  wind. 

Since  the  7-mph  wind  recorded  for  the  Johnie  Boy  event  had 
no  discernible  effect  on  the  deposition  pattern,  a  large 
part  of  the  minus  No.  200  material  found  xq)on  the  collectors 
in  the  outer  rings  must  have  been  deposited  by  the  base 
surge  or  by  the  later  effects  of  weather. 

Figure  4.2,  which  gives  the  settling  velocities  of 
quartz  in  air,  shows  that  fragments  larger  than  100  to 
200  mm  might  not  be  accelerated  enough  by  the  explosion  to 
reach  terminal  velocities  and  would  therefore  behave  bal- 
listically.  Because  of  their  high  terminal  velocities 
(caused  by  a  low  drag  force/mass  ratio),  fragments  above 
2  ram  in  diameter  could  be  expected  to  be  projected  to  much 
greater  ranges  compared  to  smaller  particles.  In  still  air, 
particles  finer  than  0.1  mm  settle  so  slowly  that  any  tiir- 
bulence  in  the  air  would  tend  to  keep  them  airborne,  and  by 
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this  means  they  could  be  transported  some  distance  from  the 
crater. 

With  the  present  limited  knowledge  on  this  subject,  it 
would  be  difficult  to  apply  the  equations  for  trajectories 
of  particles  in  air  (Appendix  B),  largely  because  of  the 
uncertainties  concerning  initial  velocities  and  air  drag. 
Also,  the  quantitative  effects  of  the  expanding  gas  bubble 
and  the  blast  wind  on  ejecta  distribution  are  not  known.  A 
final,  meaningful  evalTjation  of  ejecta  parameters  will  re¬ 
quire  utilization  of  equations  of  the  type  given  in  Appen¬ 
dix  B  and  an  extension  of  the  theory  of  trajectories  and  the 
mechanics  of  motion  of  both  large  and  small  particles  in 
air. 


4.3  DISTRIBUTION  OF  EJECTA 

The  manner  in  -which  ejecta  is  deposited  upon  the  ground 
s-urface  per  -unit  area  (5  ,  kg/m*^)  as  the  result  of  a  surface 
or  -undergro-und  explosion  is  dependent  upon:  (l)  the  yield 
or  equivalent  weight  of  explosive,  (2)  the  depth  (position) 
at  which  the  explosion  occurs,  and  (3)  the  characteristics 
of  the  mediimi  cratered  (Reference  lO).  There  are  other  fac¬ 
tors  that  influence  to  a  relatively  minor  degree  the  amo-unt 
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and  pattern  of  throwoirt,  e.g.  wind  conditions,  type  of  ex¬ 
plosive  (which  in  turn  determines  the  energy-cotipling  char¬ 
acteristics  of  the  explosion),  inhomogeneities  within  the 
medi\3ni  itself,  etc.  Until  more  data  are  available  from  a 
variety  of  nuclear  explosions  in  a  variety  of  materials, 
the  problem  may  well  be  simplified  to 

S  (R)  =  f  (W,  Z,  a)  (4.1) 

where  6  (r)  is  a  functional  relation  between  the  depo¬ 
sition  per  unit  8-rea  and  R  ,  the  range  from  GZ;  W  is  the 
yield  (weight  of  explosive);  Z  is  the  scaled  position  of 
the  charge  with  reference  to  the  air-ground  interface;  and 
o  is  the  strength  properties  of  the  medium  cratered.  For 
the  Johnie  Boy  shot,  W  ,  Z  ,  and  o  are  simply  parameters 
of  the  experiment:  hence,  the  deposition  per  unit  area  de¬ 
pends  mainly  upon  the  range  from  GZ,  or 

8  =  f  (R)  (4.2) 

The  observed  pattern  of  ejecta  deposition  for  Johnie 
Boy  is  shown  by  means  of  a  contoured  map  in  Figure  3*2.  In 
general,  the  Johnie  Boy  distribution  pattern  conforms  to 
those  obtained  from  similar  experiments  conducted  wider 
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similar  weather  conditions.  For  Johnie  Boy,  quantities  of 
ejecta  deposition  per  unit  area  large  enough  to  be  of  con- 
cern  as  an  overbiu*den  problem  (say  5  kg/m  ,  or  approximately 
1  Ib/ft  )  generally  occurred  within  a  radial  distance  of 
less  than  200  meters,  or  within  about  10  apparent  crater 
radii  of  GZ. 

In  Figure  4.3,  ejecta  distribution  for  Johnie  Boy  is 
plotted  logarithmically  as  a  function  of  radial  distance 
from  GZ.  Average  values  of  ejecta  thicknesses  were  em¬ 
ployed  with  an  assumed  specific  gravity  of  1.5  to  calculate 
the  ejecta  mass  (weight)  inside  the  50-meter  ring,  while 
the  mean  values  of  ejecta  weights  on  recoverable  collectors 
were  used  to  determine  ejecta  distribution  at  radial  dis¬ 
tances  corresponding  to  collector  rings  (Table  4.l).  The 
total  mass  of  ejecta,  ,  was  then  calculated  from  the 
eq\iation 

/oo 

6  (R)  R  dR  (4.3) 

o 

where  5  (r)  is  the  equation  of  the  curve  shown  in  Figure 
4.3  as  determined  by  the  method  of  least  sqxiares.  The 
Johnie  Boy  ejecta  data  are  thus  described  by 

6  =  2.06  X  10®  (4.4) 
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The  mass  (weight)  of  ejecta  between  the  lip  crest  and  infin¬ 
ity  was  determined  from  Equations  4.3  and  4.4  to  be  approxi¬ 
mately  4.45  X  10^  kg. 

4.4  EJECTA  EQUATIONS 

The  results  of  a  number  of  ejecta  studies  are  summa¬ 
rized  in  Appendix  A.  In  each  case,  the  amoxmt  of  ejecta 
deposited  per  imit  area  is  expressed  by  an  equation  of  the 
form  of  Eq\aation  4.4,  or 

6  =  k  r""  (4.5) 

where  k  is  a  constant  for  a  given  set  of  data,  R  is  the 
radial  distance  from  GZ  in  meters,  and  m  is  the  slope 
(characteristically  negative)  of  the  c\irve  assuming  a 
straight  line  on  logarithmic  paper. 

There  is  some  question  as  to  the  total  applicability  of 
an  equation  of  the  form  of  Equation  4.5.  For  example,  if 
the  Johnie  Boy  ejecta  data  are  plotted  on  semilogarithmic 
paper,  the  plot  in  Figure  4.4  resiilts.  This  plot  indicates 
that  there  are  two  regimes  of  exponential  behavior.  One 
regime  applies  over  the  range  from  the  peak  of  the  crater 
lip,  which  occurs  at  a  distance  of  about  23  meters  from  GZ, 
out  to  approximately  50  meters.  The  second  regime  applies 
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over  the  interval  from  about  50  meters  out  to  the  most  dis¬ 
tant  station  from  which  ejecta  was  recovered,  or  at  least 
600  meters. 

The  superiority  of  one  system  of  graphical  presenta¬ 
tion  of  the  data  over  the  other  (compare  Figures  4.3  and 
4.4)  is  not  clearly  established.  While  the  data  appear  to 
fit  a  hyperbolic  formulation  equally  well,  there  are  theo¬ 
retical  s\aggestions  that  an  exponential  formulation  would  be 
more  appropriate,  e.g.  Newton's  law  of  cooling,  which  repre¬ 
sents  q\iite  well  almost  any  phenomenon  that  gradually  dimin¬ 
ishes,  as  does  the  thickness  of  throwout  deposition. 

In  order  to  compare  Johnie  Boy  ejecta  results  with  re- 
siilts  of  other  studies,  a  hyperbolic  formulation  is  used  in 
this  report.  Futwe  work  in  this  field  should  examine  in 
more  detail  the  exponential  approach,  particularly  as  re¬ 
gards  the  two  or  more  regimes  of  behavior.  An  extensive 
effort  to  further  develop  this  approach  at  this  time  is  not 
possible  because  data  are  not  available  in  sufficient  detail. 

4.5  PREDICTION  OF  EJECTA  DEPOSITION 

The  ejecta  equations  presented  in  Table  A.l  (Appen¬ 
dix  a)  can  be  transposed  to  the  dimensionless  form 
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(Reference  lO)  and  written  as  follows: 


6  r  d 
a  a 

W 


=  b(  ^ 
a 


m 


(4.6) 


Solution  for  the  coefficient  b 

k  r  d 
_ a  a 

W 


yields  the  following: 
m 


=  b 


b  = 


,  l+m  , 
k  r  d 
a _ a 

W 


(4.7) 


Table  4.2  presents  the  transposed  ejecta  equations  (in  di¬ 
mensionless  form)  for  specific  shots  in  desert  alluvium. 
Figures  4.5  and.  4.6  provide  plots  of  these  equations.  If 
the  Stagecoach  II  curve  is  omitted,  primarily  because  of  its 
abnormal  slope,  the  characteristic  slope  shown  in  Figure  4.6 
represents  the  average  slope  of  all  remaining  ctirves.  In 
effect,  this  asstmies  that  the  characteristic  slope  for 
ejecta  curves  in  desert  alluvium  can  be  adequately  repre¬ 
sented  by  a  function 


Such 

tion 


6  r  d  /  „  \-3.4 
^  a. 

a  supposition  pres\mies  that  the  factor  b 
of  the  depth  of  the  explosion.  Figure  4.7 


(4.8) 

is  a  func- 
presents  a 
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plot  of  b  versus  scaled  depth  of  burst,  Z'  .  The  correla¬ 
tion  appears  to  be  adequate  for  both  HE  and  NE  charges. 
Stagecoach  II  again  was  omitted  from  the  plot  for  the  reason 
stated  above. 

Ejecta  deposition  can  be  thus  predicted  for  a  wide  vari¬ 
ety  of  depths  of  btirst  by  using  Figure  4.7  in  conjunction 
with  Equation  4.8  and  the  predicted  apparent  crater  radius 
and  depth. 

4.6  REIATIVE  DISTRIBUTION  OF  EJECTA. 

The  fraction  of  total  ejecta  volume  (dimensionless)  as 
a  function  of  range  from  GZ  in  units  of  crater  radii  is 
shown  in  Figure  4.8.  The  Johnie  Boy  shot  fits  well  within 
the  envelope  of  data  established  by  previous  HE  smd  NE  shots 
in  desert  aJ-luviion.  Figure  4.8  shows  that  53  percent  of  the 
Johnie  Boy  ejecta  was  deposited  within  two  crater  radii  from 
GZ,  75  percent  within  three  crater  radii,  and  approximately 
90  percent  within  five  crater  radii.  Thus^most  of  the 
ejecta  is  deposited  within  ranges  wherein  only  hardened  fa¬ 
cilities  will  be  expected  to  survive. 

An  equation  for  estimating  the  thickness  of  ejecta 
deposition  (t^)  for  various  explosion  yields  in  desert 
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alluvium  (surface  or  near- surface  bvirst)  can  be  developed  as 
follows.  Since  the  ejecta  deposition  6  (kg/m*")  can  be 
defined  as 


o  =  c  ( 

^  ■‘a 


m 


(4.9) 


and  the  ejecta  thickness  can  be  determined  by  dividing  6 
by  the  specific  weight  (density)  of  the  material,  it  follows 
that 


or 


t  =  - 
e  7 


(rj 

^  a  ^ 

The  scaled  ejecta  thickness  is  then  defined  by 

(r 

>  Q 


(4.10) 


^e  c  ^ 

'  7  W^/3.4 


_  (4.11) 

7  ^a  ^ 

Using  the  ejecta  data  for  Jangle  HE-2  and  for  Johnie  Boy 
(tabulated  with  that  for  other  desert  alluvium  shots  in 
Table  4.3),  Equation  4.11  becomes 


t  =  5-3  X  10"^  ^  A  ^  ^ 

e  \  r 

a 


(4.12) 


Equation  4.12  is  applicable  to  surface  (or  very  near 
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surface)  explosions  in  desert  alluvium  and  can  be  used  for 
estimating  the  ejecta  thickness  in  meters. 


4.7  REIATION  BETWEEN  EJECTA  VOLUME  AND  CRATER  VOLUME 

A  correlation  between  the  total  volume  of  ejecta  (V^) 
and  the  apparent  crater  volume  (V  )  for  various  shots  in 
desert  eilluvium  is  shown  in  Table  4.4  and  Figure  4.9-  The 
relation  is  satisfactorily  described  by 

V,  .  V/  (4.13) 

where  is  a  constant  and  n  an  exponent.  Since  the  ex- 

J. 

ponent  of  V  in  Equation  4.13  is  less  than  1,  the  ratio 

3. 

V  /v^  decreases  as  V  becomes  larger,  indicating  that  the 
volimie  within  the  apparent  crater  attributable  to  compres¬ 
sion  and  plastic  flow  increases  as  the  apparent  volimie 
increases  but  not  in  direct  proportion. 

If  the  volimie  due  to  gross  compression  and  plastic  flow 
is  considered  to  be  equal  to  the  apparent  crater  volume  less 
the  total  ejecta  volvme,  then 
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TABLE  4.1  AVERAGE  EJECTA  DATA  FOR  JOHNIE  BOY 


Station 

R,  Radial 
Distance 
from  GZ 

8,  Weight  of 
Ejecta  per 
Itoit  Area 

6  r 

a 

d  * 
a 

r 

a 

w 

meters 

kg/m^ 

Lip 

24.4 

1,890 

X 

CVJ 

10"^ 

1.31 

Lip 

27.4 

537 

2.06  X 

10"^ 

1.47 

Lip 

36.6 

350 

1.34  X 

10"^ 

1.97 

Lip 

45.7 

231 

8.84  X 

10-2 

2.46 

Ring  A 

50 

13.3 

5-09  X 

10"^ 

2.69 

Ring  B 

75 

10.7 

4.10  X 

10"^ 

4.03 

Ring  C 

115 

8.37 

3.22  X 

10"^ 

6.18 

Ring  D 

175 

2.52 

9.62  X 

10-*^ 

9.41 

Ring  E 

265 

0.39 

1.49  X 

10-*^ 

14.3 

Ring  F 

400 

0.11 

4.21  X 

10"^ 

21.5 

Ring  G 

600 

0.004 

1.53  X 

10-^ 

32.2 

*  r  = 
a 

18.6  meters,  d 

a 

=9*3  meters. 

and  W  =  4.53  X 

10^  kg 
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TABLE  h.3  EJECTA  EQUATIONS  FOR  EXPLOSIONS  AT  SIMILAR  SCAIED  DEPTHS  OF  BURSTS 
IN  DESERT  ALLUVIUM 


Event 

Explo¬ 

sive 

Yield,  W 

&  -  Ejecta  Distribution,  kg/m^ 

leg 

kR® 

c(R/r^)“ 

Near-Surface  Bursts 

1.71  X  10®  (R/r^^)"®'®’^ 

Jangle  HE-2 

HE 

1.81  X  10^ 

7.6  X  10^  R’^'^'^ 

Johnie  Boy 

NE 

4.53  X  10^ 

2.06  X  10® 

4.7T  X  10®  (R/r^^)”®'^^ 

Denths  of  Burst 

-  65  ft/kW^ 

6.85  X  10*^ 

Teapot  Ess 

NE 

1.1  X  10^ 

2.1  X  10^ 

Stagecoach  II 

HE 

1.81  X  10*^ 

1.95  X  10^  r"^’^ 

2.34  X  10^ 

ERA  112 

HE 

1.16  X  10^ 

9.46  X  10® 

2.40  X  10®  (R/r^)"®’”^^ 

ERA  U5 

HE 

1.81  X  10*^ 

2.05  X  lo”^  R"®'^'^ 

1.0  X  10®  (R/rg^)"®’^"^ 

Depths  of  Burst 

;  ~  150  ft/kt^/^ 

2.54  X  10*^  (r/i-^)"®’^ 

Sedan 

NE 

9.1  X  10^ 

5.77  X  10^  R‘®'®** 

Scooter 

HE 

4.5  X  10^ 

7.95  X  10®  r"®’®° 

2.46  X  10®  (R/r^)"®’®° 

Stagecoach  III 

HE 

1.81  X  10*^ 

7.32  X  lo"^  r"®'®^ 

4.24  X  10®  (R/r^)“®‘®^ 
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*  Source  material:  References  1  sind  8. 
Assumed  density  =1-5  gm/cmS 


DISTANCE  FROM  GZ,  METERS 

Figiire  4.1  Variation  of  screen  size  fractions  with  distance 
from  GZ.  Note  ring  designations  on  abscissa.  Samples  F-1 
and  G-32  are  not  considered  in  this  graph. 
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PARTICUE  SIZE,  MICRONS 


Figure  4.2  Settling  velocities  of  quartz  in  e 
(Reference  15). 
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Figtire  4.3  Ejecta  distribution,  Johnie  Boy.  Data 
points  represent  arithmetic  means^of  recoverable  sam¬ 
ples.  One  data  point,  1,890  kg/m  at  24.4  meters,  is 
not  shown.  Data  points  for  rings  A  and  B  were  not  used 
in  computation  of  curve. 
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Figure  4.8  Fraction  of  total  ejecta  volume  as  a  function  of  range  from  GZ. 
Envelope  is  from  Reference  2. 
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CHAPTER  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  CONCLUSIONS 

The  Erode  mathematical  model  (Reference  U)  appears  to 
serve  quite  well  as  a  qualitative  explanation  of  the  move¬ 
ment  of  material  from  craters  in  alluvium,  althoiigh  it  was 
derived  to  apply  to  tuff. 

A  straight-line,  logarithmic  plot  of  5  versus  R  ap¬ 
pears  to  give  a  good  representation  of  ejecta  distribution 
between  the  crest  of  the  crater  lip  and  the  outer  limit  to 
which  throwout  is  a  problem  as  such  (i.e.  excluding  fallout). 
This  treatment  also  provides  a  means  of  comparing  the  vari¬ 
ous  ejecta  studies  that  have  been  conducted  to  date. 

Within  the  range  of  yields  for  which  experience  is 
available  in  ejecta  studies,  prediction  of  throwout  is 
feasible  by  the  methods  developed  in  Chapter  4.  For  desert 
alluvium,  90  percent  or  more  of  all  ejecta  can  be  expected 
to  fall  within  10  crater  radii  of  GZ.  The  Johnie  Boy  shot, 
which  should  closely  approximate  a  low-yield,  tactical 
ground  burst,  indicates  that  throwout  should  not  be  an  over¬ 
burden  problem  to  hardened  military  installations  beyond 
this  distance. 
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Particle  size  analysis  showed  that  there  was  consider¬ 
able  sorting  of  ejecta  by  air  drag,  the  percentage  of  mate¬ 
rial  finer  than  screen  size  No.  10  decreasing  rapidly  with 
distance  from  GZ.  Much  of  the  finest  material,  especially, 
minus  No.  200,  is  believed  to  have  been  deposited  by  the 
base  surge  or  by  postshct  weather.  Wind  did  not  influence 
that  part  of  the  ejecta  which  vrould  constitute  an  overburden 
problem. 

In  Johnie  Boy  Project  1.5,  a  relatively  simple  experi¬ 
ment  provided  good  data  on  and  a  better  insight  into  the 
problem  of  nuclear  surface  bursts,  such  as  might  be  expected 
in  the  tactical  employment  of  nuclear  weapons.  The  experi¬ 
mental  technique  herein  reported,  including  both  sample  col¬ 
lection  and  close-in  excavation,  may  be  improved  by  experi¬ 
ence  and  more  adequate  funding. 

5.2  RECOMMENDATIONS 

Prediction  of  throwout  phenomena  for  extremely  large 
yields  will  require  additional  experimentation.  Also,  re¬ 
finement  of  the  knowledge  gained  thus  far  in  this  field  is 
desirable.  For  these  purposes,  full  advantage  should  be 
taken  of  large-scale  explosive  testing,  especially  shallow 
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HE  and  NE  biirsts  in  the  kt  range  and  similarly  positioned 
HE  bursts  in  the  Mt  range. 

A  better  technique  of  ejecta  sampling  is  desirable.  In 
addition  to  the  use  of  collectors  for  fine  particles  and  ex¬ 
cavation  through  the  crater  lip  to  define  close-in  ejecta, 
the  preshot  spraying  of  areas  beyond  the  predicted  crater 
lip  with  a  colored,  cohesive  paint  is  suggested.  In  con- 
Jmction  with  some  soil;  of  a  coordinate  system,  this  should 
permit  the  taking  (or  simply  measuring  in  place)  of  many 
more  ejecta  samples,  or  more  important,  of  larger  samples  at 
greater  distances,  which  would  give  a  more  accurate  repre¬ 
sentation  of  ejecta  distribution  and  total  mass. 

Futiire  ejecta  studies  should  include  a  consideration 
of  an  exponential-  formixlation,  as  discussed  in  Chapter  4, 
and  shovild,  if  possible,  define  the  regimes  characterizing 
throwout  behavior. 

Finally,  the  initiation  of  a  program  to  establish  math¬ 
ematical  models  for  trajectories  of  ejecta  with  varying  par¬ 
ticle  size  distributions,  and  taking  into  account  air  drag 
and  wind  velocities,  is  also  recommended. 
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APPENDIX  A 

SUMMARY  OF  CRATERIWG-EJECIA  EXPERIMENTS 

Table  A.l  simniarizes  all  usable  cratering-ejecta  experi¬ 
ments  known  to  have  been  conducted  to  date,  including 
Johnie  Boy. 
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tabu:  A.l  CRATER  EJECTA  STUDIES 


Shot  Detiena- 
tloo  or  Evant 

Yield 

Depth  of 
Burat 

Z* 

■>1 

■■a 

Ejecta  Distribution 
Equation 

Cratered. 

Material 

Kllotooa  (kt) 
or  Tone  (T) 

kg 

M 

./4«^* 

• 

■ 

n 

.  .3 

.* 

kg/«® 

RE 

Sedan 

100  kt 

9.1  X  10'^ 

193 

0.43 

97 

- 

195 

5.07  X  10* 

- 

5.77  X  10**  R"*'"* 

Desert  alluvlun 

Teapot  £•• 

1.2  kt 

1.1  X  10^ 

20.3 

0.20 

27. k 

- 

kk.k 

7.36  X  10* 

— 

2.1  X  10**  R**'%* 

Desert  alluvlua 

Danny  Boy 

0.ii2  kt 

1.26  X  10^ 

33.5 

0.45 

19.2 

**3-5 

32.6 

3.06  X  10** 

— 

5.69  X  10^ 

Basalt 

Johole  Boy 

0.50  kt 

1.53  X  10’ 

0.6 

0.01 

9-3 

-- 

18.6 

k.O?  X  10^ 

-- 

2.06  X  10® 

Desert  alluviun 

HE 

Scooter 

0.5  kt 

4.5  X  10* 

38 

0.50 

22.8 

.. 

k6.6  ' 

—  7.4  X  10* 

- 

7.95  X  10®  R'*'*° 

Desert  alluviun 

Stagecoach  II 

20  T 

1.8  X  10“ 

5.2 

0.20 

7.2 

-- 

15. k 

—  2.37  X  10* 

~ 

1.95  X  10®  R-*-*® 

Desert  alluviun 

Stagecoach  HI 

20  T 

1.8  X  10* 

10.4 

0.40 

8.6 

-- 

17.8 

—  4.11  X  10* 

- 

7-32  X  10^  R“^'^ 

Desert  alluviun 

SurrieLl 

5  T 

1.5  X  10* 

Contact 

-0.0 

3-8 

- 

6.1 

.. 

- 

4.5  X  10*  R"*'** 

Silt  and  clay 

Suffleld 

100  T 

9.1  X  lo'' 

Contact 

-0.0 

6.2 

- 

23.5 

—  3-55  4  10* 

- 

2.56  X  10®  R"*'®" 

Silt  and  clay 

White  Tribe  I>III 

6  T 

5.2  X  10* 

Contact 

-0.0 

2.5 

- 

5.3 

- 

8.9  X  10*  R"*'** 

Caliche 

Jangle  KE>2 

20  T 

1.8  X  lo‘‘ 

1.6 

0.05 

k.5 

- 

12.1 

—  9-91  X  10* 

- 

7.6  X  10® 

Desert  alluviun 

Buckboard  (11) 

20  T 

1.8  X  10*^ 

7.8 

0.30 

7.6 

- 

13.6 

--  1.53  X  10^ 

- 

2.42  X  10^  R"*'*® 

Basalt 

Buckboard  (12) 

20  T 

1.8  X  10*^ 

13.0 

0.50 

10.6 

- 

17. k 

--  3-82  X  10* 

- 

1.35  X  10®  R-*-®® 

Basalt 

Buckboard  (13) 

20  T 

1.8  X  lo‘‘ 

17.9 

0.68 

k.9 

- 

U.2 

6.56  X  10® 

- 

3.88  X  10*  R-*-®® 

Basalt 

ERA  603 

0.033  T 

I.I15  X  10* 

0.8 

0.15 

- 

1.5 

3-5  .  — 

2.26  X  10* 

2.21  X  10*  B**  ''* 

Sandstone 

ERA  807 

0.033  T 

l.*>5  X  10* 

0.8 

0.15 

-- 

1.6 

-- 

k  .k  •• 

4.13  X  10* 

6.23  X  10*  R'*'** 

Sandstone 

ERA  806 

0.033  T 

1.15  X  10* 

0.8 

0-15 

- 

1.8 

•• 

k.O  •• 

2.89  X  10* 

7.02  X  10*  R"*'  '* 

Sandstone 

ERA  819 

0.033  T 

1.45  X  10* 

0.8 

0.15 

- 

2.0 

- 

4.8 

k.07  X  10^ 

2.60  X  10*  R"*’*' 

Sandstone 

ERA  80U 

0.033  T 

1.45  X  10* 

1-5 

0.29 

- 

2.3 

k.3 

4.07  X  10* 

1.51  X  10*  R-*’®* 

Sandstone 

ERA  605 

0.033  T 

1.45  X  10* 

3.8 

0.72 

- 

k.5 

- 

2.8 

3.37  X  10* 

3.20  X  10*  B"*’** 

Sandstone 

ERA  610 

1.26  T 

1.16  X  10* 

1.5 

0.1k 

- 

3.0 

-- 

9.9 

2.44  X  10* 

8.83  X  10*  B"*  *''' 

Saxtdstone 

ERA  8U 

1.26  T 

l.l£  X  10* 

1.5 

0.1k 

-- 

3.2 

-- 

7.7 

1.99  X  10* 

2.T1  X  10*  B'*’'*® 

Sandstone 

ERA  SlU 

20  T 

1.8  X  10* 

3.8 

0.14 

- 

8.2 

- 

17.2 

3. 06  X  10* 

3.59  X  10®  B-*-®® 

Sandstone 

ERA  615 

20  T 

1.8  X  10* 

3.3 

0.1k 

- 

82 

-- 

21.5 

3. S'*  X  10* 

6.10  X  10*  B‘*’*® 

Sandstone 

ERA  617 

160  T 

1.45  X  10* 

7.6 

O.IO 

Ik. 3 

-- 

28.9 

1.45  X  10* 

1.51  X  10?  B"*'** 

Sandstone 

ERA  501 

0.03  T 

1.45  X  10* 

2.0 

O.ko 

- 

2.8 

- 

7.25  X  10* 

5.94  X  10*  B'*'®* 

Limestone 

ERA  UO 

0.16  T 

1.45  X  10* 

1.1 

0.21 

2.3 

2.7 

k.O 

5.1  4.53  X  10* 

7.92  X  10* 

3.83  X  10*  B-*-9® 

Dry  sand 

ERA  106 

1.28  T 

1.16  X  10^ 

0.8 

0.10 

3.0 

3-2 

5.6 

8.2  l.kT  X  10® 

3.11  X  10* 

3.6k  X  10**  R"®*^^ 

Dry  sand 

ERA  U2 

1.28  T 

1.16  X  10^ 

2.1 

0.20 

3.7 

5.3 

9.1 

10.3  3.68  X  10* 

7.06  X  10* 

9.46  X  10®  b"*'** 

Dry  sand 

ERA  115 

20  T 

li 

1.81  X  10 

0.21 

7.0 

9.1 

22.9 

25.2  5.91  X  10^ 

8. 49  X  10* 

2.05  lo'’ 

Dry  sand 

ERA  316 

0.03  T 

5.0  X  10 

0.7 

0.20 

1.8 

- 

2.7 

—  2.09  X  lO* 

- 

4.96  X  10*  B-*-*® 

Dry  clay 

ERA  30I1 

0.16  T 

l,ti5  X  10^ 

l.l 

0.20 

1.8 

3.1 

3.2 

4.2  2.32  X  10* 

8.78  X  10* 

7.45  X  10** 

Dry  clay 

ERA  310 

0.16  T 

1.45  X  10* 

1.1 

0.20 

2.1 

3.1 

3.k 

k.k  2.5k  X  10^ 

6.21  X  10* 

1.10  X  10^ 

Dry  clay 

ERA  313 

0.16  T 

1.45  X  10* 

1.1 

0.20 

2.k 

2.9 

3.9 

k.7  k.25  X  10^ 

8.50  X  10* 

2.k0  X  10^ 

Dry  clay 

ERA  305 

0.16  T 

I.U5  X  10® 

2.1 

o.ko 

2.1 

3.5 

3-6 

5.0  3.68  X  10^ 

1.30  X  10* 

k.35  X  10^  R"^'®^ 

Dry  clay 

ERA  Synaaetry 

0.16  T 

I.U5  X  10® 

2.1 

o.ko 

2.1 

3.2 

3.8 

4.6  3.68  X  10* 

9.90  X  10* 

4.02  X  10**  r“^'^ 

Dry  clay 

ERA  30B 

1.28  T 

1.16  X  10* 

0.8 

0.10 

3.7 

k.l 

6.1 

7.5  1.52  X  10® 

2.83  X  10* 

1.22  X  10* 

Dry  clay 

ERA  309 

1.28  T 

1.16  X  kP 

2.1 

o.ko 

k.7 

5.3 

6.6 

9.0  2.0k  X  10® 

5.38  X  10* 

1.12  X  10® 

Dry  clay 

ERA  312 

1.28  T 

1. 16  X  10* 

2.1 

o.ko 

k.6 

5.5 

7-9 

9.3  3.68  X  10® 

7.08  X  10* 

1.12  X  10*  B-*-** 

Dry  clay 

ERA  315 

20  T 

1.81  X  10** 

5.3 

o.ko 

12.8 

Ik. 3 

19.5 

23.6  5.38  X  10* 

9.90  X  10^ 

1.62  X  10®  B"*''*'' 

Dry  clay 

ERA  318 

160  T 

1.45  X  10^ 

10.7 

0.2 

1B.3 

23-5 

36.6 

39.6  3-12  X  10* 

u 

5.10  X  10 

9.91  X  10^ 

Dry  clay 

ERA  ‘‘02 

0.16  T 

1.45  X  10® 

O.B 

0.15 

3.0 

- 

5.7 

7.2  1.16  X  10® 

- 
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— 
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APPENDIX  B* 


PARTICLE  TRAJECTORIES 


B.l  TRAJECTORIES  IN  A  VACUUM 

The  well-knowii  parametric  equations  for  the  parabolic 
trajectory  of  a  projectile  in  a  gravitational  field,  neg¬ 
lecting  the  effects  of  the  atmosphere,  are  as  follows: 

(B.l) 
(B.2) 

Figiore  B.l  illustrates  these  equations.  Eliminating  t 


X  =  V  t  sin  6 
o  o 

4-  ^  1  4.2 

y  =  v^  t  cos  0^  -  g  t 


and  solving  the  resulting  expression  for  v^  give  v^  as 

a  function  of  x  ,  y  ,  and  0^  .  If  an  observed  pair  of 

values  X  =  R  and  y  =  h  are  inserted,  v  as  a  function 
p  ’  o 


of  0  is  given  by 
o 


1/2 


°  (R  sin  2  0  -  2  h  sin^  9 

P 


(B.3) 


v_  = 

fR  sin  20  -  2  h  sin^  (  , 

o  o' 

If  two  different  pairs  of  observed  values  of  x  and  y  are 


inserted,  two  independent  equations  which  uniquely  define 


V  and  0  result, 
o  o 


*  Prepared  by  G.  B.  Clark  and  T.  J.  Kinzer,  Jr.,  during 
their  summer  employment  at  WES,  1963. 
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B.2  DRAG  ON  PARTICLES  IK  A  FLUID 


As  soon  as  a  medium  is  introduced,  drag  forces  should 
be  considered.  The  treatment  given  this  phenomenon  is  gov¬ 
erned  by  the  physical  problem  being  considered.  Given  an 
airplane  of  known  design,  the  question  concerns  the  details 
of  airflow  around  it  at  any  velocity.  The  problem  here  is 
the  distribution  of  throwout  from  an  explosion.  Even  if  the 
exact  aerodynamic  behavior  of  some  given  particle  could  be 
predicted  (and  it  cannot),  the  details  would  be  useless  in 
view  of  the  great  variety  of  particle  shapes  which  comprise 
the  ejecta.  Thus  the  simple  approach  by  which  the  assump¬ 
tion  is  made  that  the  ejecta  may  be  characterized  by  certain 
empirical  parameters  inserted  into  the  theory  for  drag  on  a 
sphere  is  an  acceptable  step  toward  producing  some  sort  of 
distribution  f\anction. 

If  temperature  effects  are  considered  to  be  restricted 
to  density  effects,  the  drag  on  a  sphere  should  depend  upon 
its  size  and  velocity  and  on  the  density  and  viscosity  of 
the  mediiim.  Dimensional  homogeneity  then  gives 


D  =  K 


,n  2-n  n-1  n 

d  n  p 


(B.4) 


If  this  is  to  be  applied  to  arbitrary  particles,  the 
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parameters  k  and  n  must  be  kncvm.  The  parameter  k  can 
be  considered  as  depending  upon  the  particle  shape;  n  de¬ 
pends  upon  velocity,  but  fort\mately  it  usually  varies 
slowly.  Conventionally,  the  streamline  (laminar)  flow  re¬ 
gime  is  associated  with  n  =  1  and  the  turbulent  regime 
with  n  =  2'  .  These  values  will  be  used  here.  It  should  be 
noted,  however,  that  the  drag  on  a  projectile  at  high  speed 
is  not  exactly  proportional  to  v  (Reference  16).  In  the 
intermediate  regime,  the  flow  breaks  away  from  the  particle 
and  changes  considerably.  In  this  regime,  n  can  be  less 
than  one,  and  equations  using  n  =  2/3  will  be  shown  here 
to  characterize  this  region.  In  fact,  n  varies  continu¬ 
ously,  and  in  the  cases  of  very  high  Reynolds  numbers 

C 

(>  5  X  10"^ ),  there  are  velocities  at  which  a  small  increase 
in  V  actually  reduces  drag. 

The  form  of  the  equations  of  motion  for  each  of  these 
values  of  n  will  be  given.  Their  derivation  follows. 

Let  the  position  of  a  particle  be  specified  by 

r‘=Tx  +  j'y+kz  (B.5) 

in  the  coordinate  system  shown  in  Figure  B.2.  Let  its 
velocity  with  respect  to  the  ground  be 

u=r'=ix  +  jy  +  kz  (B.6) 
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Then  its  acceleration  with  respect  to  the  grotind  is 

u  =  ^  =  "i  X  +  y  +  k  z  (B.7) 

Assume  that  a  constant  horizontal  wind  of  velocity 

blows  at  an  angle  a  with  respect  to  the  z-axis,  i.e.  in 

the  ^  =  a  direction.  The  components  of  the  wind  are 

V  =  (V  sin  cc)  i  +  (v  cos  a)  k  (B.8) 

w'w  '  'w 

..'k 

The  drag  force  is  determined  by  the  particle  velocity  v 
with  respect  to  the  air,  however,  not  the  ground,  and  this 
velocity  is 


V  =  u  -  V 


w 


Cb.9) 


Since  drag  is  opposite  in  direction  from  particle  velocity, 
the  drag  force  can  be  written  as 


2-n 

n-1 

1 

-K 

d 

Po 

1  V 

2-n 

n-1 

-K 

d 

P 

1  V 

in  V 

Ivl 
in-l  - 


(B.IO) 


Gravitational  force  is  directed  vertically  downward,  and  for 
a  particle  of  density  p  in  a  medium  of  density  ,  it  is 
given  by 

Fg  =  (p  -  p^)  g  (B.ll) 
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3 

where  a^d  is  the  voluiae  of  the  particle.  This 
written  as 

F  =  -M*  g 

O 

where  M'  is  the  effective  mass  of  the  particle, 
ancy  is  negligible,  then  M'  =  M  .  Since  the  net 
the  s\mi  of  the  component  forces, 


M  r'  =  D  +  ? 

g 


=  -K 


2-n 


o 


-  M’ 


g  J 


Examination  of  Figure  B.2  shows  that 
X  =  r  sin  6  sin  ^ 
y  =  r  cos  6 
z  =  r  sin  9  cos  ^ 

Thus, 


V  =  (u  sin  6  sin  ^  sin  a)  i  +  u  cos 


+  (u  sin  9  cos  0  -  V  cos  cc)  k 

and 


2  --  -- 
V  -  V  •  V 


=  u^  -  2  u  V  sin  9  cos  -  a)  +  V  ^ 
w  ^  w 
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(B.12) 
If  buoy- 
force  is 


(B.13) 

(B.lU) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 
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Further, 

.  (v=)~ 


r 

=  -  2  u  V  sin  6  cos  (^  -  a)  +  V  ^ 

1  W 


Cb.19) 


and  finally, 


Ti  „  jH  2-n  n-1  2  r.  IT  •  n  fd 

r=-/fdii  p  u-2uV  sin  0  cos  (0  -  a) 

^  o  w  / 


n-1 

+  ^  sin  0  sin  ^  sin  o:)  i 

+  u  cos  0  j  +  (u  sin  0  cos  ^  cos  a)  if] 


-  M'  g  J 


(B.20) 


The  corresponding  orthogonal  components  can  now  he  separated 
to  give  three  simultaneous  equations.  For  n  =  2  , 

M  X  =  -K  d^  j^u^  -  2  u  sin  0  cos  -  a) 

+  V  ^  sin  0  sin  ^  -  V  sin  oc)  (B.21) 

w  J  w 

M  y  =  -/c  d^  p^  ^u^  -  2  u  sin  0  cos  {0  -  a.) 

+  V  ^  (u  cos  0)  -  M’  g  (B.22) 

w 
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2  2 

M  z  =  -/c  d  pQ  ^  ~  ^  ^  ® 

+  (u  sin  6  cos  ^  -  y^  cos  a) 

For  n  =  1  , 

M  X  =  -/<  d  u  (u  sin  6  sin  0  -  V  sin  a) 

'  V 

My  =  -/<d^^-M•  g 
M  z  =  -/f  d  |ji  (u  sin  Q  cos  0-1  cos  a) 

'  V 


For  n  =  2/3  , 

=  -K  (1^/^  ^  ^  ^  ~ 


M  X 


+  V 


V 


■1/6 


(u  sin  e  sin  0-1^  sin  a) 


M  y  =  -K  ■  2  u  V  sin  9  cos  {0  - 


U  'i  =  -K  d' 


2/3  „V3  „  -1/3 


u 


2  u  V  sin  9  cos  (/  - 


+  V  ^  fV6  Q  j,Qg  0  .  Y  COB  a) 

w  J  ^  w 

Making  use  of  the  relations  in  Equations  B.l4  thro\agh 
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(B.23) 

(B.24) 

(B.25) 

(B.26) 

a) 

(B.27) 

a) 

(B.28) 

a) 

(B.29) 


4 


rnMcincKiTi  A  I 


(B.37) 


(b.38) 

If  the  wind  velocity.  V  ,  is  zero  and  the  motion  is  put  in 

w 

the  xy -plane,  these  sets  of  equations  become; 


For 

n  = 

2  , 

M  X 

=  -K 

2 

^  Po 

.2 

X 

sin 

e 

(B.39) 

M  y 

=  -K 

d"  P 

.2 

_JL_ 

cos 

0 

-  M’ 

g 

(b.4o) 

For 

n  = 

1  , 

M  X 

=  -K 

d  p  X 

(b.4i) 

M  y 

=  -K 

d  h  y 

-  M' 

g 

(B.42) 

For 

n  = 

2/3  , 

M  X 

=  -K 

4/3 

Po 

-1/3 

(sin  0)^/3 

(b.43) 

M  y 

=  -K 

,V3 

Po 

-1/3 

(cos  0)^/2  y^/2  _  M'  g 

(b.44) 

M  z 

in  all  of 

these 

cases  is  zero. 

M  y  =  ,V3  ^^-V3 


r  .2 


Lcos^  6 


-  2  y  tan  6  cos  (0  -  a)  +  y  -  M'  g 


M  z  =  -K  d^/^  p 


--V3  U 


.2 

z 


2  2 
sin  e  cos  0 


2  z  V 


cos 


^  cos  (^  -  a)  +  (z  -  cos  a) 
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In  the  case  n  =  1  ,  the  Equations  B.33  throijgh  B.35 
and  B.41  through  B.42  can  be  integrated;  however,  this  is  a 
special  case,  and  for  other  values  of  n  the  equations  are 
coupled  and  nonlinear.  Therefore,  approximate  methods  are 
reqiiired.  For  small,  heavy  particles  the  n  =  2  regime  is 
the  most  important.  For  the  condition  of  no  wind  and  verti¬ 
cal  motion  only.  Equations  B.30  through  B.32  or  B.39  through 
B.40  can  be  integrated  in  closed  form.  For  two  approximate 
solutions,  see  Section  B.l|. 

B.3  TERMINAL  VELOCITIES 

The  relation  M  r  =  D  +  F  implies  that  when  D  =  F 

s  s 

the  acceleration  is  zero,  leaving  velocity  constant.  This 
equlllbrl\am  velocity  is  called  the  terminal  velocity.  The 
vector  nature  of  the  defining  relation  means  that  terminal 
velocity  is  confined  to  vertical  motion;  however,  the  nxmier- 
ical  value  representing  the  terminal  velocity  of  a  given 
particle  in  a  given  medium  is  useful  in  cases  of  nonvertical 
motion . 

For  vertical  motion,  the  force  due  to  gravity  is 

5^  =  (p  -  Pq)  g  J  (B.U) 
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and  the  drag  force  is 


n  2-n  n-1  |->-|n-l  -*  /t, 

D  =  -K  d  1  v|  V  (B.IO) 

The  terminal  velocities  for  the  three  values  of  n  con¬ 
sidered  in  this  appendix  are  as  follows: 

For  n  =  2  (tiurhulent ), 


For  n  =  1  (streamline). 


For  n  =  2/3  (intermediate). 


=  K, 


(B.45) 


(B.46) 


(B.47) 


The  k's  are  dependent  upon  particle  shape.  The  effective 
diameter  d^  in  B.47  Is  empirical  and  is  given  hy 

d^  =  d  -  ^  d  (B.48) 


where  d  is  the  true  diameter,  d  is  given  hy 


d  =' 


V  36.^  , 

e  Po  (p  -  Pq) 


(B.49) 


and  ^  has  been  erapiricainy  determined  to  be  0.4  for 
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spheres  and  0.28  for  irreg\ilar  quartz  particles,  expressed 
in  centimeter-gram-second  units. 

E:q)erimental  vork  on  terminal  velocities  of  irregular 
crushed  quartz  particles  (Reference  17 )  has  shown  that  par¬ 
ticles  smaller  than  about  25  microns  obey  the  streamline 
flow  conditions,  while  those  between  85  and  2,000  microns 
exhibit  terminal  velocities  which  vary  linearly  with  diame¬ 
ter.  Larger  particles  show  turbulent  flow.  For  these 
larger  particles,  resistance  increases  proportional  to  the 
sq\iare  of  the  diameter  and  weight  increases  proportional  to 
the  cube  of  the  diameter i  consequently,  drag  becomes  rela¬ 
tively  less  important. 

Very  fine  particles  of  dleuneter  smaller  than  the  mean 
free  path  of  air  molecules  require  a  correction  to  the  ob¬ 
served  velocity  to  calculate  drag.  This  correction  is 


where  K  is  a  constant  (0.86  for  air)  and  X  is  the  mean 
free  path  of  the  air  molecules. 


b.4  approximate  solutions  to  the  equations  of  motion 


Since  the  case  n  =  1  is  the  only  one  for  which  the 
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equations  of  motion  can  be  integrated  in  closed  form,  ap¬ 
proximate  procedures  are  necessary  for  the  more  normal  case 


i 


in  which  this  value  of  n  does  not  apply.  Typically, 
these  involve  numerical  methods,  assumptions  to  lineetrize 
or  decouple  the  eqiiations,  or  both.  The  results  are  then 
tabulated  in  parametric  form,  as,  for  example,  in  artillery 
firing  tables. 

For  firing  tables,  one  method  used  is  the  following: 

Let  be  the  angle  between  the  particle's  velocity 

vector  and  the  horizontal.  Then 
d0 

cos  ’  ®1  degrees  (B.51) 


This  approximate  relation  and  the  assumption  that  the  di*ag 
coefficient  is  a  constant  give 


X  =  K  (  sec^  )  0^ 


(B.52) 


or 


dv  d0 

=  X  sec^  0  -rrr- 
dt  1  dt 


which  yields 

V 


(f  "  t)}]„  °  (B-53) 
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This  can  be  solved  for  numerical  values  by  the  method  of 
q,Tiadrat\ares . 

Another  approximation  for  entirely  different  purposes 
is  given  in  Reference  l8.  Suppose  a  vertically  falling 
raindrop  at  terminal  velocity  falls  into  a  layer  of  air 
with  a  horizontal  motion  perpendiciilar  to  the  raindrop's 
path.  The  assumption  that  the  raindrop's  vertical  motion 
remains  unchanged  gives  the  following  results: 

Let  v^  =  vertical  velocity  =  constant 

v^  =  horizontal  velocity  of  raindrop 

D  =  «:  v^  (B.54) 

dv  / 


Integration  of  Equation  B.55  from  initial  velocity  v^^  to 


v^  yields 

t 


(B.56) 


where  t  is  the  elapsed  time  between  the  two  limits. 
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dt 


to  =  't  (-4-^  ^ 


exp  2 


(B.5T) 


•where 


A  =  In 


''xl  -I 


(B.58) 


The  horizontal  distance  traveled  in  time  t  is  foimd  from 
a  further  integration  to  be 


(B.59) 
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Figure  B.l  Coordinate  system  for  vacuum  trajectory, 


Figure  B.2  Coordinate  system  for  general 
trajectory  analysis. 
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APPENDIX  C* 


GEOLOGIC  OBSERVATIONS  OF  THE  JOHNIE  BOY  CRATER  LIP 

C.l  GENERAL 

Three  trenches,  one  each  on  the  north,  south,  and  east 
sides  of  the  crater  (see  Figure  2.8  of  main  text),  were  ex¬ 
cavated  down  to  or  below  original  groiand  surface.  These 
trenches,  trending  along  radials  projected  from  GZ,  extended 
to  the  approximate  limit  of  the  crater  lip. 

C.2  CRITERIA  FOR  RECOGNITION  OF  ORIGINAL  GROUND  SURFACE 
Observations  of  the  walls  of  the  trenches  indicated 
that,  althotigh  there  was  no  pronounced  stratification  of  ma¬ 
terials,  there  were  slight  and  subtle  differences  in  lithol¬ 
ogy,  color,  and  consistency  which  permitted  the  recognition 
of  several  "beds"  or  "layers"  within  the  alluvium.  These 
layers,  present  in  all  trenches,  were  numbered  in  consecu¬ 
tive  order  from  uppermost  (youngest)  to  lowermost  (oldest) 
and  are  illustrated  in  Figures  C.l  throu^  C.5»  The  atti¬ 
tude  and  sequence  of  these  layers  were  used  to  determine  the 
amount  of  uplift  and  deformation  caused  by  the  explosion. 

*  Prepared  by  R.  W.  Hunt,  R.  T.  Saucier,  and  W.  B. 
Steinrlede,  Jr.,  WES  Soils  Division. 
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The  general  location  of  the  post shot  (uplifted  origi¬ 
nal)  grovind  surface  vas  determined  hy  Inference  from  the 
sequence  of  layers  of  alluvitun.  The  precise  location  of 
this  surface  could  be  obtained  at  various  points  in  the 
trench  walls  by  noting  the  locations  of  grass,  roots,  stems, 
and  occasionally  fragments  of  cloth  or  other  types  of  litter. 

C.3  TRMCH  OBSERVATIONS 

C.3.1  East  Wall  of  North  Trench.  As  shown  in  Figiire 
C.l,  there  were  two  sequences  of  beds  or  layers  separated 
by  a  shear  plane  which  dipped  away  from  the  apparent  crater 
at  an  angle  of  about  40  degrees.  The  postshot  ground  sur¬ 
face  had  been  displaced  upward  at  least  4  feet  near  the  lip 
crest  and  had  deformed  in  the  vicinity  of  the  shear  plane. 

The  warping  of  the  postshot  groimd  surface  just  north 
of  the  shear  plane  and  the  configuration  of  the  layers  in 
the  vicinity  of  the  shear  plane  suggest  that  the  northern¬ 
most  sequence  was  uplifted  and  tightly  folded  before  the 
shear  plane  developed.  The  shear  plane  appears  to  be  the 
result  of  the  innermost  ( southernmost )  sequence  being 
thrust  outward  and  beneath  the  outermost  sequence. 

C.3»2  East  Wall  of  South  Trench.  In  the  approximately 
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l6  to  l8  ft  of  trench  wall  not  obscured  by  talus,  it  was 
possible  to  define  a  situation  in  which  the  alluvial  se¬ 
quence  had  been  completely  overturned  or  folded  over  (Figure 
C.2).  The  three  layers  of  alluvium  which  occurred  below  the 
original  ground  surface  were  present  in  inverted  order  over- 
lying  a  zone  in  which  the  order  was  still  preserved  but 
slightly  uplifted.  The  fold  axis  was  preserved  and  is  ap¬ 
parent  both  in  the  schematic  diagram  and  in  the  accompanying 
photograph.  The  material  between  the  two  segments  of  origi¬ 
nal  ground  surface  was  a  mixture  of  loose,  sandy  silt, 
gravel,  and  grass  roots,  and  it  probably  represented  a 
deposition  of  ejecta  which  occurred  before  the  foldover 
was  completed. 

C.3-3  West  Wall  of  South  Trench.  A  small,  talus-fr?e 
exposure  in  this  wall  revealed  that  the  alluvial  sequence 
had  been  uplifted  to  an  angle  of  about  65  to  70  degrees  (see 
Figure  C.3).  The  beds  occurred  about  4  to  5  feet  higher 
than  they  did  in  the  opposite  wall  of  the  trench,  which  was 
about  15  feet  to  the  east,  and  they  actually  protruded  about 
6  inches  above  the  ejecta  surface. 

Because  of  the  obscuring  ejecta  and  talus  cover,  it  was 
not  possible  to  determine  if  foldover  occurred  at  this 
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point.  There  was  no  indication  of  a  repeated  inverted  se¬ 
quence  of  layers  lying  beyond  the  mapped  sequence  (in  a 
direction  away  from  GZ);  the  high  angle  to  which  the  layers 
were  iq)lifted  does  stiggest,  however,  that  foldover  may  have 
occiorred . 

C.3.4  West  Wall  of  North  Trench.  A  small  exposure  in 
this  trench  wall  (Figure  C.U)  exhibited  a  normal  sequence  of 
alluvial  layers  overlain  by  an  inverted  sequence .  The  two 
segments  of  the  original  ground  surface  were  in  contact, 
with  no  separating  layer  of  ejecta  as  shown  in  Figure  C.2. 

In  this  instance,  the  repeated  sequence  of  beds  is 
strong--but  not  conclusive — evidence  that  foldover  occurred. 
Since  an  actual  fold  axis  was  not  identified,  there  is  a 
possibility  that  the  inverted  sequence  represented  a  large 
block  of  alluvium  which  was  ejected  during  the  explosion 
and  landed  intact  but  inverted. 

C.3.5  Spilth  Wall  of  East  Trench.  An  arciiate  exposure 
here  revealed  two  steeply  tilted  sequences  of  alluvium:  an 
innermost  sequence  in  normal  order  and  an  outermost  sequence 
in  inverted  order  (Figure  C.5)-  The  two  segments  of  origi¬ 
nal  ground  surface  were  in  direct  contact  with  one  another 
and  dipped  away  from  the  crater  at  approximately  70  degrees. 
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The  lip  crest  of  this  was  formed  by  a  mass  of  layer  2 

which  protruded  through  the  ejecta. 

Althoiagh  the  absence  of  a  fold  axis  again  prevented  the 
positive  identification  of  foldover,  the  high  angle  uo  which 
the  beds  or  layers  were  tilted  and  the  length  along  which 
the  back-to-back  ground  siorfaces  could  be  traced  (about 
7  feet)  are  strong  indications  that  foldover  occurred. 

C.4  GENERAL  OBSERVATIONS  AROUND  CRATER 

The  upper  slopes  of  the  apparent  crater  were  covered 
with  fallback.  This  material,  consisting  of  loose,  tan 
silt,  sand,  gravel,  and  cobbles,  masked  for  the  most  part 
the  uplifted  beds  or  layers  of  alluvium.  Around  almost  the 
entire  crater,  the  lip  crest  was  characterized  by  steep, 
outward-dipping  exposures  of  layered  alluvium.  The  general 
appearance  was  that  of  an  irregtilar,  roughly  circular 
"hogback"  ridge.  At  each  of  the  measured  sections,  the  only 
layers  exposed  in  this  ridge  were  layers  1  and  2.  In  a1 T 
cases,  layer  3  was  covered  with  ejecta  or  talus.  The  rather 
steeply  dipping  uppermost  10  to  15  feet  of  the  ejecta  sur¬ 
face  (immediately  outward  from  the  lip  crest)  probably  di¬ 
rectly  reflected  the  underlying  tilted  alluvial  sequences. 
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The  angle  of  repose  of  the  ejecta  in  this  area  was 
occasionally  as  high  as  35  to  40  degrees. 


C.5  CONCLUSIONS 

Observations  in  the  trenches  and  around  the  crater  in¬ 
dicate  that  in  all  cases  the  originally  horizontaJ.  or  nearly 
horizontal  layers  of  alluvium  were  displaced  upward  and  out¬ 
ward  by  the  force  of  the  explosion.  Although  warping  of 
beds  and  layers  may  have  been  most  common,  in  at  least  one 
case  (Figure  C.l)  a  large  block  of  alluvium  was  thrust  be¬ 
neath  another  along  a  well-defined  shear  plane . 

In  addition  to  the  uplifting,  or  upthrust ing,  of 
layers,  there  were  several  examples  of  complete  foldover, 
or  overturning  of  layers.  Additional  excavation  around 
the  crater  would  probably  reveal  that  this  foldover  was  a 
rather  common  occxirrence. 

Because  of  the  discontinuous  occurrences  of  foldover 
and  the  shearing  and  distortion  of  the  alluvium,  it  is  not 
possible  to  accurately  predict  the  behavior  and  degree  of 
upthrust  of  the  original  ground  surface  which  will  occur 
as  the  result  of  a  near- surface  detonation  such  as  the 
Johnie  Boy  event. 
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Figure  C.l  Photograph  and  schematic  diagram  of  exposed 
beds  in  east  vidll  of  north  trench. 
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Figure  C.2  Photograph  and  schematic 
beds  in  east  wall  of  south  trench. 
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APPROXIMATE  SCALE  IN  FEET 
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Figure  C.3  Photograph  and  schematic  diagram  of  exposed 
beds  in  west  wall  of  south  trench. 
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Figure  C.4  Photograph  and  schematic 
beds  in  west  wall  of  north  trench. 
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APPROXIMATE  SCALE  IN  PEET 
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APPROXIMATE  SCALE  IN  PEET  PROM  ASSUMED  POINT 


Figure  C.5  Schematic  diagram  of  exposed  beds  in  south  wall 
of  east  trench. 
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